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November 7, 2024 at 12:32 


1. Intro. This program was derived directly from XCCDC, by adding proof logging. We will use the 
VeriPB program as the external verifier. This program creates two extra output files: an OPB file containing 
a system of linear inequalities to model our XCC problem, and a PBP file which contains our proof logging. 
These proof logging extensions were added by Filip Stappers in 2024. 

This program is an experimental XCC solver, which often looks ahead considerably further than DLX2 
does. More precisely, it maintains “domain consistency”: An option O is eliminated when its use would cause 
some primary item I ¢ O to have no options remaining. In a sense, I’m performing the work of DLX-PRE 
repeatedly as the search proceeds. With luck, the total time will decrease, although the time per node is 
potentially much larger. 

Furthermore, I’m continuing to experiment with sparse-set data structures, as I did in the similar program 
SSXCC, which was inspired by Christine Solnon’s XCC-WITH-DANCING-CELLS. 

The DLX input format used in previous solvers is adopted here, without change, so that fair comparisons 
can be made. (See the program DLX2 for definitions. Much of the code from that program is used to parse 
the input for this one.) 
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2. After this program finds all solutions, it normally prints their total number on stderr, together with 
statistics about how many nodes were in the search tree, and how many “updates” were made. The running 
time in “mems” is also reported, together with the approximate number of bytes needed for data storage. 
(An “update” is the removal of an option from its item list, or the removal of a satisfied color constraint 
from its option. One “mem” essentially means a memory access to a 64-bit word. The reported totals don’t 
include the time or space needed to parse the input or to format the output.) 

Empirical tests show that this program takes more elapsed time per mem than most other programs that 
I’ve written. I don’t know why. Perhaps it’s because the number of “global registers” is unusually large. 

Here is the overall structure: 


7##define o memst++ /* count one mem */ 

#define 00 mems += 2 /* count two mems */ 

##define 000 mems +=3 /* count three mems */ 

#tdefine subroutine.overhead mems += 4 

#define O "%" /* used for percent signs in format strings «/ 
#define mod % /* used for percent signs denoting remainder in C */ 
7##define maz_stage 500 /* at most this many options in a solution */ 


7##define maz_level 50000 /* at most this many levels in the search tree */ 
#define maz_cols 10000 /* at most this many items «/ 


7#define maz_nodes 50000000 /* at most this many nonzero elements in the matrix */ 
#define poolsize 100000000 /*x at most this many entries in pool */ 
##define savesize 1000000 /* at most this many entries on savestack */ 


#-define bufsize (9 * maz-_cols + 3) /* a buffer big enough to hold all item names «/ 


#include <stdio.h> 

#include <stdlib.h> 

#include <string.h> 

#include <ctype.h> 
typedef unsigned int uint; /* a convenient abbreviation */ 
typedef unsigned long long ullng; /* ditto */ 


(Type definitions 9); 

(Global variables 4); 

(Subroutines 7); 

int main(int argc,char *argv[]){ register int c, cc, i, 3, k, p, pp, q, r, 8, t, cur-choice, best_itm, 

line , optnr; 
(Process the command line 5); 
(Input the item names 29); 
(Input the options 31); 
(Output our problem to the OPB file for proof logging 40); 
(Output header of the PBP file for proof logging 45); 
if (vbose & show_basics) (Report the successful completion of the input phase 38 ); 
if (vbose & show_tots) (Report the item totals 39); 
imems = mems,mems = 0; 
if (baditem) (Report an uncoverable item 37) 
else (Solve the problem 59); 
done: (Say adieu 6); 


} 
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You can control the amount of output, as well as certain properties of the algorithm, by specifying 


options on the command line: 


‘vy (integer )’ enables or disables various kinds of verbose output on stderr, given by binary codes such as 
show_choices; 


e ‘m( integer )’ causes every mth solution to be output (the default is mO, which merely counts them); 
e ‘d( integer )’ sets delta, which causes periodic state reports on stderr after the algorithm has performed 


approximately delta mems since the previous report (default 10000000000); 


e ‘c( positive integer )’ limits the levels on which choices are shown during verbose tracing; 
e ‘C( positive integer )’ limits the levels on which choices are shown in the periodic state reports; 


‘1(nonnegative integer )’ gives a lower limit, relative to the maximum level so far achieved, to the levels 
on which choices are shown during verbose tracing; 

‘t (positive integer )’ causes the program to stop after this many solutions have been found; 

(integer )’ sets timeout (which causes abrupt termination if mems > timeout at the beginning of a level); 
(filename )’ (required) name of the “.opb file” and “.pbp file” for proof logging; 

(filename )’ to output a “shape file” that encodes the search tree; 

x( positive integer )’ causes partial solutions of this many stages to be written to files, not actually explored; 
‘X( filename )’ to input and resume a partial solution. 


‘T 
“Pp 
‘S 


‘ 


7##define show_basics 1 /* vubose code for basic stats; this is the default */ 

#define show-_choices 2 /* vbose code for backtrack logging */ 

#define show_details 4 /* vubose code for stats about choices «/ 

#define show_purges 8 /* vbose code to show inconsistent options deleted «/ 
##define show_supports 16 /* vbose code to show new supports */ 

##define show_option_counts 32 /* vbose code to count active options */ 

#define show_mstats 64 /* vbose code to show memory usage in key arrays */ 
#define show_profile 128 /*x vbose code to show the search tree profile «/ 

#define show_fullstate 256 /* vbose code for complete state reports */ 

#define show_tots 512 /* vbose code for reporting item totals at start */ 

#define show_warnings 1024 /* vbose code for reporting options without primaries */ 
#t#define show_maz_deg 2048 /* vbose code for reporting maximum branching degree */ 
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4. (Global variables 4) = 
int vbose = show_basics + show_warnings; /* level of verbosity */ 
int spacing; /* solution k is output if k is a multiple of spacing */ 
int show_choices.maz = 1000000; /* above this level, show_choices is ignored */ 
int show_choices.gap = 1000000; /* below level mazl — show-choices.gap, show_details is ignored */ 
int show_levels.max = 1000000; /* above this level, state reports stop */ 
int mazl, mazs; /* maximum level and stage actually reached «/ 
int xcutoff = —1, xcount; /* stage when partial solutions output, and their number «/ 
int mazsaveptr ; /* maximum size of savestack */ 
char buf [bufsize]; /* input buffer «/ 
ullng count; /* solutions found so far */ 
ullng options; /* options seen so far */ 
ullng imems, mems, bmems, nmems, pmems, tmems; /* mem counts */ 
ullng updates; /* update counts */ 
ullng bytes; /* memory used by main data structures */ 
ullng nodes; /* total number of branch nodes initiated */ 
ullng thresh = 10000000000; /* report when mems exceeds this, if delta 40 */ 
ullng delta = 10000000000; /* report every delta or so mems */ 
ullng maxcount = *ffffffffffffffFfFf; /* stop after finding this many solutions «/ 
ullng temeout = *1fffffffffffffff; /* give up after this many mems */ 
FILE xshape_file ; /* file for optional output of search tree shape */ 


char *shape_name; /* its name */ 

FILE xopb_file; /* file for OPB model of the problem we are solving */ 
char opb_name [255]; /* its name */ 

FILE «pbp_file; /* file for our proof log */ 

char pbp_name [255]; /* its name «/ 

int maxdeg; /* the largest branching degree seen so far */ 

int plLconstraints; /*x Current number of proof logging constraints written */ 


See also sections 10, 19, 76, and 96. 


This code is used in section 2. 
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5. If an option appears more than once on the command line, the first appearance takes precedence. 
(Process the command line 5) = 
for (j = argc —1,k =0; 7; j——) 
switch (argv[j][0]) { 


case ’v’: k |= (sscanf (arguv|j] + 1,""O"d", &ubose) — 1); break; 

case ’m’: k |= (sscanf (arguv|j] +.1,""O"d", &spacing ) — 1); break; 

case ’d’: k |= (sscanf (arguv[j] + 1,""O"11d", &delta) — 1), thresh = delta; break; 
case ’c’: k |= (sscanf (arguv[j] + 1,""O"d", &show_choices.maz ) — 1); break; 
case ’C’: k |= (sscanf (arguv[j] + 1,""O"d", &show-_levels.max ) — 1); break; 

case ?1’: k |= (sscanf (arguv[j] + 1,""O"d", &show_choices_gap ) — 1); break; 
case ’t’: k |= (sscanf (argu[j] + 1,""O"11da", &maxcount) — 1); break; 

case ’T’: k |= (sscanf (argu[j] + 1,""O"11d", &timeout) — 1); break; 

case ’S’: shape._name = argu|j] + 1, shape_file = fopen (shape_name, "w"); 


if (ashape_file ) 
forintf (stderr, "Sorry, ,Ijcan’ tyopen file,‘"O"s’ for writing! \n", shape.name); 


break; 


case ’P’: strcpy(opb_name, argu|j] + 1), strcat(opb_name,".opb"), opb_file = fopen(opb_name, "w"); 
if (sopb_file) fprintf(stderr, "Sorry, ,Ican’tyopen file, ‘"O"s’ for writing! \n", opbname); 


strepy (pbp_name, argu |j] + 1), strcat (pbp_name,".pbp"), pbp_file = fopen(pbp_name, "w"); 


if (=pbp_file) fprintf (stderr, "Sorry, ,Ican’ tyopen file, ‘"O"s’ for writing! \n", pbp_name); 


break; 
case ’x’: k |= (sscanf (arguv|j] + 1,""O"d", &xcutoff) — 1); break; 
case ’X’: (Open «cutoff_file for reading, and break 93); 
default: k = 1; /* unrecognized command-line option */ 
} 
if (k V sopb_file V spbp_file) { 
fprintf (stderr, "Usage: "O"sP<bar>, [v<n>] |, [m<n>] , [d<n>] ,[c<n>] ,[C<n>] " 
"(1<n>] ,[t<n>] ,[T<n>] , [S<bar>] ,[x<n>] ,[X<bar>] |<foo.d1x\n", argu [0]); 
exit (—1); 
} 


This code is used in section 2. 
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6. Idon’t report the memory used for deg, stagelevel, and profile, because they are only for documentation, 
not part of the search process. 
(Say adieu 6) = 
( Write contradiction line for proof logging 49); 
if (vbose & show_profile) (Print the profile 90); 
if (vbose & show_maa_deg) fprintf (stderr, "The, maximum branching, degree, was,,"O"d.\n", mazrdeq); 
if (vbose & show_basics) { 
forintf (stderr, "Altogether,."O"1llu,solution"O"s,"O"llu+"O"1lu,mems,", count, 
count =1?"":"s",imems,mems); 
bytes = (itemlength + setlength) * sizeof (int) + last_node * sizeof 
(node) + (4 * mars + mazl) « sizeof (int) + mazsaveptr * sizeof (twoints) + poolptr * sizeof 
(twoints); 
fprintf (stderr, "\"O"1llujupdates, "O"1llu,bytes,"O"llu,nodes,", updates, bytes, nodes ); 
fprintf (stderr, "acosty"O"11d%%, ubcosty"O"11d%%, wccosty"O"11d%%.\n", 
mems ? (200 * nmems + mems)/(2 * mems) : 0,mems ? (200 * pmems + mems)/(2 * mems) : 0, 
mems ? (200 * bmems + mems)/(2* mems) : 0); 


if (vbose & show_mstats) { 
fprintf (stderr, "\yitemlength="O0"d, ,setlength="O"d, jlast_node="O"d;\n", itemlength, setlength, 
last_node); 
forintf (stderr, "\maxsaveptr="O0"d, ,poolptr="O"d, maxstage="O"d,,maxlevel="O"d.\n", 
maxsaveptr , poolptr, maxs, mazl); 


if (sanity_checking) fprintf (stderr, "sanity_checking,was,on! \n"); 
if (leak_checking) fprintf(stderr, "leak_checking,was,jon! \n"); 

if (shape_file) fclose (shape_file ); 

if (opb_file) fclose (opb_file); 

if (pbp_file) fclose (pbp_file); 

if (acount) (Report the number of partial solutions output 95); 


This code is used in section 2. 


, ere’s a subroutine for use in debuggin u ope it’s never invoked. 
7. Here’ broutine f debugging, but I hope it’ ked 
(Subroutines 7) = 
void confusion(char id) 
{ /* an assertion has failed */ 
rintf (stderr, "troubleafter,, umems ,., unodes: js! \n", mems, nodes, id); 
printf (stderr, "troubleyafter,,"O"lld "O"1lldnodes:,%s!\n", , nodes, id); 
See also sections 12, 13, 14, 15, 16, 20, 21, 22, 23, 24, 25, 41, 50, 51, 52, 57, 58, 80, 87, 88, and 89. 


This code is used in section 2. 
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8. Data structures. Sparse-set data structures were introduced by Preston Briggs and Linda Torczon 
[ACM Letters on Programming Languages and Systems 2 (1993), 59-69], who realized that exercise 2.12 
in Aho, Hopcroft, and Ullman’s classic text The Design and Analysis of Computer Algorithms (Addison— 
Wesley, 1974) was much more than just a slick trick to avoid initializing an array. (Indeed, TAOCP exercise 
2.2.6—24 calls it the “sparse array trick.” ) 

The basic idea is amazingly simple, when specialized to the situations that we need to deal with: We can 
represent a subset S of the universe U = {2,%1,...,%n—1} by maintaining two n-element arrays p and gq, 
each of which is a permutation of {0,1,...,n — 1}, together with an integer s in the range 0 <s <n. In 
fact, p is the inverse of g; and s is the number of elements of $. The current value of the set S' is then simply 
{Lpo,+-+,Lp,_,}. (Notice that every s-element subset can be represented in s!(n — s)! ways.) 

It’s easy to test if x, € S, because that’s true if and only if gq, < s. It’s easy to insert a new element x, 
into S: Swap indices so that p, = k, q, = s, then increase s by 1. It’s easy to delete an element x, that 
belongs to S: Decrease s by 1, then swap indices so that p, = k and q, = s. And so on. 

Briggs and Torczon were interested in applications where s begins at zero and tends to remain small. In 
such cases, p and q need not be permutations: The values of ps, ps41, .--,; Pn—1 can be garbage, and the 
values of gz need be defined only when xz, € S. (Such situations correspond to Aho, Hopcroft, and Ullman, 
who started with an array full of garbage and used a sparse-set structure to remember the set of nongarbage 
cells.) Our applications are different: Each set begins equal to its intended universe, and gradually shrinks. 
In such cases, we might as well maintain inverse permutations. The basic operations go faster when we know 
in advance that we aren’t inserting an element that’s already present (nor deleting an element that isn’t). 

Many variations are possible. For example, p could be a permutation of {xo,21,...,%n—1} instead of a 
permutation of {0,1,...,2—1}. The arrays that play the role of qg in the following routines don’t have 
indices that are consecutive; they live inside of other structures. 


8 DATA STRUCTURES XCCDC-PROOF 89 


9. This program has an array called item, with one entry for each item. The value of item[k] is an 
index x into a much larger array called set. The set of all options that involve the kth item appears in that 
array beginning at set[{x]; and it continues for s consecutive entries, where s = size(a) is an abbreviation 
for set|a — 1]. If ttem[k] = x, we maintain the relation pos(x) = k, where pos(x) is an abbreviation for 
set[x — 2]. Thus item plays the role of array p, in a sparse-set data structure for the set of all currently 
active items; and pos plays the role of q. 

Suppose the Ath item x currently appears in s options. Those options are indices into nd, which is an 
array of “nodes.” Each node has four fields: itm, loc, clr, and atra. Ifa <q<a+s, let y = set[q]. This 
is essentially a pointer to a node, and we have nd[y].itm = x, nd[y].loc = q. In other words, the sequential 
list of s elements that begins at « = item|k] in the set array is the sparse-set representation of the currently 
active options that contain the kth item. The clr field nd[y].clr contains x’s color for this option. The itm 
and clr fields remain constant, once we’ve initialized everything, but the loc fields will change. The ztra 
field has special uses as we maintain domain consistency, as explained later. 

The given options are stored sequentially in the nd array, with one node per item, separated by “spacer” 
nodes. If y is the spacer node following an option with ¢ items, we have nd[y].itm = —t. If y is the spacer 
node preceding an option with t items, we have nd[y].loc = t. 

This probably sounds confusing, until you can see some code. Meanwhile, let’s take note of the invariant 
relations that hold whenever k, g, x, and y have appropriate values: 


pos(item|k]) =k; nd[set[q]].loc =q; item[pos(x)] =x; set|nd[y].loc] = y. 


(These are the analogs of the invariant relations p[q[k]] = q[p[k]] = & in the simple sparse-set scheme that 
we started with.) 

The set array contains also the item names, as well as two fields mark(x) and match(x) that are used for 
compatibility checking. (The match field is present only in secondary items.) 

We count one mem for a simultaneous access to the itm and loc fields of a node, also one for simultaneous 
access to clr and ztra. 


#define size(x) set[(x) — 1] /* number of active options of the kth item, 7 */ 
#define pos(x) set[(a) — 2] /* where that item is found in the item array */ 
#define Iname(x) set|(x) — 4] /* the first four bytes of ’s name */ 

#define rname(x)  set|(x) — 3] /* the last four bytes of x’s name */ 


#define mark(x) set|(x) — 5] /* a stamp for incompatible items «/ 

#-define match(x) set|(x) — 6] /* a required color in compatibility tests */ 
#tdefine primextra 5 /* this many extra entries of set for each primary item */ 
7##define secondeztra 6 /* and this many for each secondary item */ 

7#define mazextra 6 /* maximum of primeztra and secondextra */ 


(Type definitions 9) = 
typedef struct node_struct { 


int itm; /* the item x corresponding to this node «/ 
int loc; /* where this node resides in x’s active set */ 
int clr; /* color associated with item x in this option, if any */ 
int ztra; /* used for special purposes (see below) */ 
} node; 


See also section 11. 


This code is used in section 2. 
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10. (Global variables 4) += 


node nd[maz_nodes]; /* the master list of nodes */ 

int last_node; /* the first node in nd that’s not yet used «/ 

int item|maz_cols]; /* the master list of items */ 

int second = maz_cols; /* boundary between primary and secondary items */ 
int last_itm; /* items seen so far during input, plus 1 */ 

int set[maz_nodes + 6 * maz-_cols]; /* sets of active options for active items */ 
int itemlength; /*x number of elements used in item */ 

int setlength; /*x number of elements used in set */ 

int active; /* current number of active items */ 

int oactive; /* value of active before swapping out current-choice items */ 
int totopts; /* current number of active options */ 

int baditem; /* an item with no options, plus 1 */ 

int osecond; /* setting of second just after initial input «/ 


11. We're going to store string data (an item’s name) in the midst of the integer array set. So we’ve got 
to do some type coercion using low-level C-ness. 
(Type definitions 9) += 

typedef struct { 


int /, 7; 

} twoints; 

typedef union { 
unsigned char str|8]; /* eight one-byte characters */ 
twoints Ir; /* two four-byte integers «/ 


} stringbuf; 
stringbuf namebuf ; 


12. (Subroutines 7) += 
void printitem_name(int k, FILE «stream) 
{ 
namebuf .lr.l = Iname(k), namebuf .Ir.r = rname(k); 
forintf (stream, ""O" .8s8", namebuf .str); 


} 
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13. An option is identified by the names of the items it contains. For proof logging we use the line number 
where the option appears in the input as the identifier of the option. Here is a routine that prints an option, 
given a pointer to any of its nodes. If showid = 1, it also prints the value of opt — 1, which should be the 
location of the spacer just preceding opt. Otherwise it optionally prints the position of the option in its item 
list. 
(Subroutines 7) += 
void print_option(int opt, FILE *stream,int showpos,int showid) 
{ 
register int k, q, 2; 
x = ndlopt].itm; 
if (opt > lastnode Vx <0) { 
forintf (stderr, "Illegaloption,"O"d!\n", opt); 
return; 
} 
if (showid) fprintf (stream, ""O"d,‘", opt — 1); 
for (q = opt; ; ) { 
print_item_name(a, stream); 
if (2 > second A nd{(q].clr) fprintf (stream,":"O"c", nd[q].clr); 
qt; 
x = ndlq).itm; 
if (x <0) { 
q+= &,; 
fprintf (stream, "\,(Line,"O"d)", linenr (q)); 
x = nd[q].itm; 
} 


if (q = opt) break; 


k = nd|{q}.loc; 
if (showid) fprintf (stream, "’"); 
if (showpos > 0) fprintf (stream, "C"O" dof "O"d) \n",k — x +1, size(x)); 
else if (showpos =0) fprintf (stream,"\n"); 
} 
void prow(int p) 
{ 
print_option(p, stderr, 1,0); 
} 
void propt(int opt) 
{ /*x opt should be the spacer just before an option */ 
if (ndlopt].itm > 0) fprintf(stderr, ""O"d isn’ tan joptionid! \n", opt); 
else print_option(opt + 1, stderr ,0, 1); 


} 


§14 


14. 
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The print_option routine has a sort of inverse, which reads from buf what purports to be the description 


of an option and verifies it. 


(Subroutines 7) += 
int read_option (void) 


{ 


} 


15. 


register int k, gq, x, 3, opt; 
for (opt =0,k =1; 0, buf [k] > °0’ A buf [k] < °9?; k++) opt =10* opt + buf [k] — 70’; 
if ((0, buf [k] A ’u’) V (0, buf [kK +1] 4 ’°?) V (0, buf [k+ 2] A ’u’)) return —1; 
for (k += 3,q = opt +1; 0, (a = nd[q].itm) > 0; q+) { 
oo, namebuf .lr.l = Iname (x), namebuf .lr.r = rname (x); 
for (j = 0; 7 <8; j++) { 
if (anamebuf .str[j]) break; 
if (0, namebuf .str[j] # buf [k + j]) return —1; 


} 

k+= 9; /* we’ve verified the item name */ 

if (0, nd[q].clr) { 
if ((0, buf [k] A ?:’) V (0, (unsigned char) buf [k + 1] 4 nd[q].cir)) return —1; 
k += 2; 


if (0, buf [k++] 4 ’’) return —1; 


} 
if (buf [k] A ’\’’) return —1,; 
return opt +1; 


When I’m debugging, I might want to look at one of the current item lists. 


(Subroutines 7) += 
void printitm (int c) 


{ 


register int p; 
if (c < primextra V c > setlength V pos(c) < 0V pos(c) > itemlength V item|pos(c)| 4c) { 
forintf (stderr, "Illegalyitem,"O"d!\n",c); 
return; 
} 
forintf (stderr, "Item"); 
print.item_name(c, stderr); 
if (c < second) fprintf (stderr, "yC"O"dyofy"O"d) , ylength,"O"d:\n", pos(c) + 1, active, size(c)); 
else if (pos(c) > active) 
fprintf (stderr, "\,(secondary,"O"d, jpurified) ,,length,"O"d:\n", pos(c) + 1, size(c)); 
else fprintf(stderr, "\(secondary,"O"d) , ,length,"O"d:\n", pos(c) + 1, size(c)); 
for (p=c; p<c+ size(c); p++) prow(set[p)); 
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16. Speaking of debugging, here’s a routine to check if redundant parts of our data structure have gone 
awry. 
7##define sanity_checking 0 /* set this to 1 if you suspect a bug */ 


(Subroutines 7) += 
void sanity(void) 


register int k, x, 7, 1, r, q, 94; 
for (k =0; k < itemlength; k++) { 
x = item|k]; 
if (pos (i) #k) { 
fprintf (stderr, "Bad. posyfield of item"); 
print.item_name (a, stderr); 


fprintf (stderr, "\("O"d,"O"d) !\n",k, x); 


for (¢ =0; i < last_node; i++) { 
1 = ndli].itm,r = nd{i].loc; 
if (1<0) { 
if (nd[it+r+1]).itm ¢ —r) fprintf(stderr, "Bad spacerin,nodes,"O"d, "O"d!\n",i,i+r+1); 
qq = 0; 
else { 
if (l>r) fprintf(stderr, "itm>loc in node,"O"d!\n", 4); 
else { 
if (set[r] £2) { 
fprintf (stderr, "Bad locyfield for option,,"O"d of yitem",r —1+ 1); 
print_item_name (1, stderr); 
fprintf (stderr, "\yinynode,"O"d! \n", 2); 


a 


if (pos(l) < active) { 
if (r <1+ size(l)) q=+1; else ¢=-1; /* in or out? */ 
if (q* qq <0) { 
forintf (stderr, "Flipped status_atoption,,"O"d of item",r —1+ 1); 
print_item_name (I, stderr); 
forintf (stderr, "\yinnode,,"O"d!\n", 2); 


qq = | 
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17. Domain consistency. The data structures above were fine for SSXCC, but this program aims to 
prune its search tree by maintaining “domain consistency.” Several more things are therefore needed. 

We regard the given XCC problem as a special case of the general binary CSP, in which the variables 
are the primary items. The domain of primary item p is the set of options that contain p. And there’s a 
constraint between each pair of primary items p and p’: Option o for p is allowed together with option o’ 
for p’ if and only if o and o’ are compatible, in the sense that they’re either equal or they have no items in 
common, except for secondary items with identical nonzero colors. 

What does domain consistency mean in this context? “For every p 4 p’ and every o in the domain of p, 
there’s a compatible option o’ in the domain of p’.” Stating this another way, suppose o is an option. Then 
the action of choosing o, in order to “cover” its primary items, must not “wipe out” the domain of any 
primary item that’s not in o. 

When an option doesn’t meet this criterion, we remove it from consideration, thus simplifying the problem. 
The removal of an option also makes other options potentially removable. Eventually we either remove the 
last option from some item’s domain, in which case there’s no solution, or we reach a stable situation where 
all domains are nonempty and consistent. In the latter case, we’ll choose an option, for an item that has 
comparatively few of them, and we'll recursively explore the consequences of either using that option or not. 

To maintain domain consistency we shall combine the ideas of Christian Bessiére’s algorithm AC-6 
[Artificial Intelligence 65 (1994), 179-190] with Christophe Lecoutre and Fred Hemery’s algorithm AC3rm 
(IJCAI Proceedings 20 (2007), 125-130], by maintaining a table of supports: This program will essentially 
construct an array S[o,p], with an element for every option o and every primary item p, such that S/o, p] 
is equal to o’ for some compatible option o’ such that p € o’, whenever p ¢ 0; and S[o, p] = # when p € o. 
This array provides witnesses to the fact that the current domains are indeed consistent. 


18. We don’t, however, actually represent the support array S directly. Instead, we represent the inverse 
function: For each option o’, we maintain a list of all the pairs (0, p) such that S[o, p] = o’. This list is called 
the trigger list of o', because we use it to maintain the support conditions: If option o’ is removed for any 
reason, thereby leaving one or more holes in the S array, the removal will trigger a series of events that will 
refill those holes, one by one. 

Each option o also has a fizit list, containing all pairs (o', p) for which the event (0, p) has been triggered 
by o’ but the corresponding hole hasn’t yet been refilled. 

There’s also a queue Q, containing all the options o for which at least one hole currently exists. 

All of these lists — the triggers, the fixits, and the queue — are singly linked, in a array called pool, whose 
elements have two fields called info and link in familiar fashion. The trigger lists and fixit lists are stacks 
(last-in-first-out); the queue is first-in-first-out. 
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19. Internally, an option o is represented by the index of the spacer just preceding that option in nd. An 
item 2, whether primary or secondary, is represented by the index where the main data for J appears in set. 
A link is represented by its index in pool. 

Variables gfront and qrear are the indices of the front and rear of Q. More precisely, gfront points to the 
front element, the node that will be removed first; grear points to a “blank” node that follows the element 
that will be removed last. The queue is empty if and only if gfront = qrear. The contents of info(qrear ) 
and link (qrear) are both irrelevant; they will be filled in when a new element is enqueued and a new blank 
element is appended. 

Fortunately there’s room enough in the existing data structures of program SSXCC to store the two pointers 
that we need for each option: The top of o’s trigger stack, called trigger(o), is kept in location nd[o].clr; 
and the top of o’s fixit stack, called fizit(o), is kept in nd[o].ctra. We have fixit(o) = 0 if and only if o is not 
in the queue. 

We'll see later than every inactive option has an age, indicating when it was purged from the current 
partial solution. This value, age(o) appears in nd[o + 1].xtra. 

For proof logging, it is conventient to remember the line number of each option in the input. This value, 
linenr(o) appears in nd[o].cir of all primary items. 

(Kludge note: With these conventions, all of an option’s dynamic data has been squeezed into the four 
otherwise unused fields nd[o].clr, nd[o].ctra, and nd[o + 1].2tra. 


#define info(p) pool[p].l 
#define link(p) pool|p].r 
#define trigger(opt) nd[opt].cir /* beginning of the trigger stack */ 
#define fizit(opt) nd[opt|.ztra /* beginning of the fixit stack */ 
#define age(opt) nd[(opt) + 1].xtra /* when was this option last purged? */ 
#define linenr(opt) ndlopt].cir /* line where this option appeared in the input */ 
(Global variables 4) += 
twoints pool [poolsize]; /* where the linked lists live «/ 
int poolptr = 1; /* the first unused cell of pool */ 
int qfront, qrear; /* the front and rear of Q */ 
unsigned int curstamp; /* the current “time stamp” */ 
unsigned int biggeststamp; /* the largest time stamp used so far */ 
unsigned int compatstamp; /* another stamp, used for compatibility tests */ 
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20. A few basic primitive routines undergird all of our list processing. 
(When counting mems here, we consider avail and poolptr to be in global registers. The compiler could 
inline this code, so I don’t count any overhead for these subroutine calls.) 


#define avail pool(0].r /* head of the stack of available cells */ 


(Subroutines 7) += 
int getavail (void) 
{ /* return a pointer to an unused cell */ 
register int p; 
p = avail; 
if (p) { 
o, avail = link (p); 
return p; /* info(p) might be anything */ 
i 
if (poolptr ++ > poolsize) { 
forintf (stderr, "Pooljoverflow,,(poolsize="O"d) !\n", poolsize ); 
exit (—7); 
} 
return poolptr — 1; 


} 


void putavail (int p) 

{ /* free the single cell p */ 
o, link (p) = avail; 
avail = p; 


} 


21. Entries of a trigger list are pairs (0,p), with the cell that mentions option o linking to the cell that 
mentions primary item p. 
(Subroutines 7) += 
void print_trigger (int opt) 
{ 
register int p, q; 
fprintf (stderr, "trigger stack,foroption,,"); 
print_option(opt + 1, stderr ,0, 1); 
for (p = trigger (opt); p; p= link(q)) { 
q = link (p); 
forintf (stderr, "\"); 
if (info(p) > 0) { 
print_option(info(p) + 1, stderr, —1, 1); 
fprintf (stderr ,","); 
p = link(p); 
print_item_name (info (p), stderr); 
} else (Print a trigger hint 71); 
forintf (stderr, "\n"); 
} 
} 
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22. (Subroutines 7) += 
void print_triggers (int all) 
{ 
register int opt, jj, optp; 
for (opt =0; opt < last_node; opt += nd[opt].loc +1) { 
if (all) { 
qj = ndlopt + 1].ttm; /* 9j is opt’s first item */ 
if (nd[opt + 1].loc > 7j + size(gj)) continue; /* is opt in 4j’s set? */ 
} 
print_trigger (opt); 
} 
3 


23. Entries of a fixit list are pairs (0,p), with the cell that mentions option o linking to the cell that 
mentions primary item p. 
(Subroutines 7) += 
void print_fizit(int opt) 
{ 
register int p, q; 
forintf (stderr, "fixity ystack_foroption,"); 
print_option(opt + 1, stderr, —1,1); 
forintf (stderr ,":"); 
for (p = fixit(opt); p; p= link(q)) { 
q = link(p); 
forintf (stderr, "\,"); 
print_item_name (info (q), stderr); 
forintf (stderr, "["O"d]", info(p)); 
} 
forintf (stderr, "\n"); 


} 


24. (Subroutines 7) += 
void print_queue (void) 


register int p; 


for (p = qfront; p 4 qrear; p= link(p)) print_option(info(p) + 1, stderr, 0,1); 
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25. Linked lists are wonderful; but a single weak link can cause a catastrophic error. Therefore, when 
debugging, I want to be extra sure that this program doesn’t make any silly errors when it uses pool 
pointers. 

Furthermore, since I’m doing my own garbage collection, I want to avoid any “memory leaks” that would 
occur when I’ve forgotten to recycle a no-longer-used entry of the pool. 

The list_check routine laboriously goes through everything and makes sure that every cell less than poolptr 
currently has one and only one use. 

Warning: Do not call list_check at a busy time during which lists are being manipulated. Wait for a quiet 
time when all of the lists are supposedly stable and well-formed. 


##define leak_checking 0 /* set this nonzero if you suspect linked-list bugs */ 
#define signbit #8000000 
#define vet_and_set (I) 
{ if (2) < OV (1) 2 poolptr) { 
fprintf (stderr, "Badlink,,"O"d!\n", 1); 
return; 


if (link (1) & signbit) { 
forintf (stderr, "Double,link,,"O"d,,"O"1ld!\n", 1, mems); 
return; 
} 
link (1) B= signbit; 
} 
(Subroutines 7) += 
void list_check (int count) 
i 
register int p, t, opt; 
for (t = 0,p = avail; p; t++,p = signbit © link(p)) vet_and_set(p); 
if (count) fprintf(stderr, "avail size,"O"d\n", t); 
for (opt =0; opt < last.node; opt += nd[opt].loc +1) { 
for (p = trigger (opt); p; p = signbit ® link(p)) vet_and_set (p); 
for (p = fizit(opt); p; p = signbit @ link(p)) vet_and_set(p); 


for (p = qfront; ; p = signbit @ link(p)) { 
vet_and_set (p); 
if (p = qrear) break; 


for (p= 1; p< poolptr; p++) { 
if (link(p) & signbit) link(p) B= signbit; 
else fprintf(stderr, "Lost.jcell,"O"d!\n", p); 
} 
i 
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26. One of our main activities is to find options O’ that are compatible with a given option O. We do this 
by marking each item I of O with compatstamp, and also recording I’s color if J is secondary. Then, given 
a candidate O’, we can easily spot incompatibility. 

That idea works when I € O is equivalent to mark(I) = compatstamp. So we start with all mark fields 
equal to zero; and we increase compatstamp by 1 whenever starting this process with a new O. 

But there’s a hitch: If this testing is done 2°? times, compatstamp will “wrap around” to zero, and our 
test might be invalid. In such a case we can still guarantee success if we take the trouble to zero out all the 
mark fields again. 


##define badstamp 0 /* set this to 3, say, when initially debugging */ 
(Bump compatstamp 26) = 
if (++ compatstamp = badstamp) { 
for (ii =0; ti < itemlength; ++) 00, mark(item[ii]) = 0; 
compatstamp = 1; 
} 


This code is used in section 27. 


27. (Prepare the mark fields for testing compatibility with opt 27) = 
(Bump compatstamp 26); 
for (nn = opt +1; 0, (4 = nd[nn].itm) > 0; nn) { 
0, mark (ti) = compatstamp; 
if (4 > second) { 
if (0, nd[nn].clr) 0, match(it) = nd[nn].clr; 
else 0, match(i#) = —1, /* this won’t match any color «/ 
} 
} 


This code is used in sections 57 and 58. 


28. At the beginning of this section, nd[optp].itm is an item ii in the middle of some option O’. If O’ is 
compatible with opt, we want to reset optp so that nd[optp] is the spacer preceding O’. We use the fact 
that i isn’t present in opt. 
(If optp is compatible with opt, break 28) = 
for (qq = optp,nn = qq +1; nn 4 qq; nn++) { 
if (0, Gj = nd[nn].itm) <0) optp = nn + 4j —1, nn = optp; /* nn is a spacer */ 
else if (0, mark (jj) = compatstamp) { /* watch out, 7j is in opt */ 
if (jj < second V (0, nd[nn].clr = 0) V (0, nd[nn].clr 4 match (jj))) break; /* incompatible */ 


if (nn = qq) break; /* not incompatible */ 


This code is used in sections 57 and 58. 
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29. Inputting the matrix. Brute force is the rule in this part of the code, whose goal is to parse and 
store the input data and to check its validity. 

We use only four entries of set per item while reading the item-name line. 
#define panic(m) 

{ fprintf(stderr, ""O"s!\n"O"d:4"O".99s\n",m, p, buf); exit (—666); } 

(Input the item names 29) = 

line = 0; while (1) { 

if (fgets (buf , bufsize, stdin)) break; 

line ++ ; 
if (0, buf [p = strlen(buf) — 1] 4 ?\n’) panic("Inputylineway,jtoo, long"); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =? |? V abuf [p]) continue; /* bypass comment or blank line «*/ 
last.itm = 1; 
break; } 
if (slast_itm) panic("No,items"); 
for (; 0, buf[p]; ) { 

o, namebuf .lr.l = namebuf .lr.r = 0; 

for (j = 0; 7 < 8A (0, visspace (buf [p + j])); j++) { 

if (buf [p+ j] =?’:? V buf[p+ 9] =’ 1?) panic("Illegalcharacter,,in,item,name"); 
o, namebuf .str|j] = buf [p + J]; 


if (j =8A 7isspace (buf [p+ j])) panic("Item mame, too, long"); 
oo, Iname (lastitm < 2) = namebuf .lr.l, rname (last_itm < 2) = namebuf .lr.r; 
(Check for duplicate item name 30); 
last_itm-++; 
if (last_itm > maz_cols) panic("Too,many,items"); 
for (p += j +1; 0, isspace (buf [p]); p++) ; 
if (buf [pl =? 1?) { 
if (second 4 maz_cols) panic("Itemname,line,contains,,| twice"); 
second = last_itm, 
for (p++; 0, isspace (buf [p]); p++) ; 


if (second = last_itm) second = maz_cols; /* no secondaries actually named */ 


This code is used in section 2. 


30. (Check for duplicate item name 30) = 
for (k = lastitm —1; k; k-——) { 
if (0, Iname(k « 2) 4 namebuf.ir.l) continue; 
if (rname(k < 2) = namebuf.ir.r) break; 


if (k) panic("Duplicate,item,name" ); 


This code is used in section 29. 
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31. (Input the options 31) = 
while (1) { 
if (fgets (buf , bufsize, stdin)) break; 
line + ; 
if (0, buf [p = strlen(buf) — 1] 4 ?\n’) panic("Option line, too, long" ); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =’ |? V abuf [p]) continue; /* bypass comment or blank line */ 
i = last_node; /*x remember the spacer at the left of this option */ 
for (pp = 0; buf[p]; ) { 
o, namebuf .Ir.l = namebuf .lr.r = 0; 
for (j =0; j < 8A (0, 7isspace (buf [p + j])) A buf [p+ j] 42:73 G++) 0, namebuf .str|j] = buf [p + J]; 
if (47) panic("Empty item name"); 
if (j =8A 7isspace (buf [p+ j]) A buf [p+ 7] 4 ?:’) panic("Item mame, too,,long"); 
(Create a node for the item named in buf [p] 32); 
if (buf [p+ i)=?2°) { 
if (k > second) { 
if ((o, isspace (buf [p + 7 + 1])) V (0, aisspace (buf [p + 7 + 2]))) 
panic("Color,must bea single,jcharacter"); 
o, nd|last_node + (pp ? 0: 1)].clr = (unsigned char) buf [p+ 7 + 1]; 
pt+= 2; 
} else panic("Primary,item must jbe,uncolored"); 


for (p += 7 +1; 0, isspace(buf [p]); p++) ; 


} 
if (spp) { 
if (vbose & show_warnings) fprintf (stderr, "Option ignored,,(no.primary,jitems) :,"O"s", buf); 
while (last_node >i) { 
(Remove last_node from its item list 33); 


last_node ——; 
} else { 
o, nd[i].loc = last_node — 1; /* complete the previous spacer */ 
last_node++; /* create the next spacer */ 
if (last_node = maz_nodes) panic("Too many nodes"); 
options ++; 


o, nd[last_node].itm = i+ 1 — last_node; 


} (Initialize item 34); 
(Expand set 35); 
(Adjust nd 36); 


This code is used in section 2. 
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32. We temporarily use pos to recognize duplicate items in an option. 

This program shifts the items of an option, if necessary, so that the very first item is always primary. In 
other words, secondary items that precede the first primary item are actually stored in nd[last_node + 1]. 
(Create a node for the item named in buf [p] 32) = 

for (k = (last_itm — 1) «2; k; k -=4) { 

if (0, Iname(k) A namebuf .lr.l) continue; 
if (rname(k) = namebuf .ir.r) break; 


if (4k) panic("Unknown,item, name"); 
if (0, pos(k) > 2) panic("Duplicate item mamejin, this,option"); 


last_node ++; 
if (last_node + 1 > maaz_nodes) panic("Too many_nodes"); 
0, t = size(k); /* how many previous options have used this item? */ 


if (pp) { /* no primary items seen yet */ 

if ((k >> 2) < second) 00, pp =1,nd[i+ 1].itm = k > 2, nd[i + 1].loc = t, nd[i + 1].clr = line ; 

else 00, nd([last_node + 1].itm = k >> 2, nd[last_node + 1].loc = t, nd[last_node + 1].clr = 0; 

} else 00, nd[last_node].itm = k >> 2, nd[last_node].loc = t, nd[last_node].clr = (k >> 2) < second ? line 
: 0; 

o, size(k) =t +1, pos(k) = last_node; 


This code is used in section 31. 


33. (Remove last_node from its item list 33) = 
0, k = nd[last_node + 1].itm « 2; 
00, size(k)——, pos(k) =i-1, 


This code is used in section 31. 


34. (Initialize item 34) = 
active = itemlength = last_itm — 1; 
for (k = 0,7 = primeztra; k < itemlength; k++) 
oo, item|k] = j,j += (k+2 < second ? primeztra : secondextra) + size((k +1) « 2); 
setlength = 7 — 4; /* a decent upper bound */ 
if (second = maz-cols) osecond = active, second = j; /* no secondary items */ 
else osecond = second — 1; 


This code is used in section 31. 


35. Going from high to low, we now move the item names and sizes to their final positions (leaving room 
for the pointers into nb). 


(Expand set 35) = 

for (; k; k—) { 
0,j = item|k — 1]; 
if (k = second) second = J; /* second is now an index into set */ 
00, size(j) = size(k < 2); 
if (size(j) =OAk < osecond) baditem = k; 
0, pos(j) =k — 1; 
00, rname(j) = rname(k < 2), name(j) = Iname(k « 2); 
o, mark (j) = 0; 

} 


This code is used in section 31. 
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36. (Adjust nd 36) = 
for (k = 1; k < last_node; k++) { 


if (0, nd[k].itm <0) continue; /* skip over a spacer */ 
0,j = item|nd|k].itm — 1]; 
i= 7 + nd[k].loc; /* no mem charged because we just read nd|k].itm */ 


o, nd|k].ttm = j, nd[k].loc = 4; 
o, set |i] = k; 
} 


This code is used in section 31. 


37. (Report an uncoverable item 37) = 
{ 
if (vbose & show_choices) { 
forintf (stderr, "Item"); 
print_item_name (item|baditem — 1], stderr); 
forintf (stderr, "\jhasyno,options! \n"); 


} 


This code is used in section 2. 


38. The “number of entries” includes spacers (because DLX2 includes spacers in its reports). If you want 
to know the sum of the option lengths, just subtract the number of options. 


(Report the successful completion of the input phase 38) = 
forintf (stderr, "("O"11ld options, ,"O"d+"O"d items, .,"O"d entries,successfullyread) \n", 
options , osecond, itemlength — osecond, last_node); 


This code is used in section 2. 


39. The item lengths after input are shown (on request). But there’s little use trying to show them after 
the process is done, since they are restored somewhat blindly. (Failures of the linked-list implementation in 
DLX2 could sometimes be detected by showing the final lengths; but that reasoning no longer applies.) 


(Report the item totals 39) = 
{ 

fprintf (stderr, "Item ,totals:"); 

for (k =0; k < itemlength; k++) { 
if (k = second) fprintf(stderr ,"|"); 
fprintf (stderr, "\\"O"d", size (item [k])); 

} 

fprintf (stderr, "\n"); 


} 


This code is used in section 2. 
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40. Proof logging. We will convert the XCC problem to a pseudo-boolean problem, so it can be verified 
by VeriPB. In the pseudo-boolean problem, all variables are booleans (0 or 1). We create one variable for 
each option (named Ojineny), and one variable for each item in each option. Primary items are named 
Olinenr-variablename, secondary items also include their color. 

Then we add constraints for each option/item pair, so that items are selected if and only if the option is 
selected. We add the constraint that for each primary item, the sum of all its corresponding ’item-in-option’ 
variables is 1, to make sure every primary item is covered. Finally, we add a constraint for each pair of 
incompatible options due to secondary items. 

We keep count of the number of constraints that we output during proof logging. Linear equalities count 
double, since they are converted to two linear inequalities by VeriPB. 

No mems are counted during proof logging. 

For example, let’s consider again the example from DLX2: 


| A simple example of color controls 
ABC|XY 

A B X:0 Y:0 

AC X:1 Y:1 

X:0 Y:1 

B X:1 

C Y:1 


The first option will be transformed in these 4 lines, to ensure the option and its variables are both either 
true or false: 


1 o38_A 1 ~03 = 1; 
1 0o38_B 1 ~03 = 1; 
1 038_X_0 1 ~03 = 1; 
1 o8_Y_0 1 ~03 = 1; 


To ensure item A will be covered, this line is added: 
1 o8_A 1 04_A = 1; 


Finally, to ensure item X isn’t given two different colors by the options of lines 3 and 4, we add this 
constraint: 
1 ~o3_X_0 1 ~04X_1 >= 1; 


(Output our problem to the OPB file for proof logging 40) = 
plconstraints = 0; 
(Output constraints for options and items 42); 
(Output constraints to ensure each primary item is covered 43); 
(Output constraints to ensure secondary items are assigned only one color 44); 
if (opb_file) fclose(opb_file); 


This code is used in section 2. 


41. (Subroutines 7) += 
void print_pLitemname(FILE «stream, int optnr,int k,int clr) 
{ 
namebuf .lr.l = Iname(k), namebuf .lr.r = rname(k); 
forintf (stream, "0o"O"a_"O" .8s", optnr, namebuf .str); 
if (k > second) { 
if (clr) fprintf (stream,"_"O"c", clr); 
else fprintf (stream,"_[]"); 
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42. (Output constraints for options and items 42) = 
for (t = 0; i < last_node; i++) { 


k = nd[i].itm; 
if (k <0) optnr = linenr(i + 1); 
else { 


forintf (opb_file, "1u"); 
print_plitemname (opb-_file, optnr , k, nd|[i].clr ); 
forintf (opb_file, "W170" O"dy= 41; \n", optnr, namebuf .str, optnr ); 
pl_constraints += 2; 
t 
} 


This code is used in section 40. 


43. (Output constraints to ensure each primary item is covered 43) = 
{ 
int wz, ss, opt; 
for (k =0; k < active; k++) { 
ii = item|k]; 
if (ii > second) continue; 
for (s = ti, ss = s+ size(w); s < ss; st+) { 
opt = set|s]; 
forintf (opb_file, "14"); 
print_pLitemname ( opb_file, linenr (opt), item|k], 0); 
printf (opb_file, "."); 


forint (opb_file, "=,1;\n"); 
pl_constraints += 2; 
} 
i 


This code is used in section 40. 
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44. For the secondary items, we add a condition in the OPB file for each pair of incompatible options. 


(Output constraints to ensure secondary items are assigned only one color 44) = 
{ 
int w, s1, s2, ss, optl, opt2, opta; 
for (k =0; k < active; k++) { 
ii = item|k]; 
if (ii < second) continue; 
for (si = ti,ss = s1 + size(ti); s1 < ss; s1++) { 
opt1 = set|s1]; 
for (s2 =s1 +1; s2 < ss; s2++) { 
opt2 = set[s2]; 
if (nd[opt1].clr 4 nd[opt2].clr V nd[opt1].clr =0V nd[opt2].clr =0) { 
fprintf (opb_file,"147"); 
for (optz = opt1; nd[opta|.itm > 0; opta——) ; /* go to spacer for linenr */ 
print_pLitemname ( opb_file, linenr (optz + 1), ttem[k], nd[opt1].clr); 
fprintf (opb_file, "41u~"); 
for (optz = opt2; nd[opta|.itm > 0; opta——) ; /* go to spacer for linenr */ 
print_pLitemname ( opb_file, linenr (optz + 1), item[k], nd[opt2].clr); 
fprintf (opb_file, "W>=41;\n"); 
plconstraints + ; 


This code is used in section 40. 


45. We also need to write the necessary steps to the OPB file, so that VeriPB can verify our proof. This 
uses the routines below. 
(Output header of the PBP file for proof logging 45) = 

forintf (pbp_file, "pseudo-Boolean, proof,,version,,1.0\n"); 

forint (pbp_file, "£"O"d\n", plconstraints ); 


This code is used in section 2. 


46. (Write a solution for proof logging 46) = 
Sprintf (pbp_file, "vu"); 
int opt; 
for (k =0; k < stage; k++) { 
for (opt = choice [levelstage[k]]; nd[opt].itm > 0; opt—) ; 
optnr = linenr (++ opt); 
forintf (pbp_file, "0" O"d,", optnr); 
for (opt; nd[opt].itm > 0; opt++) { 
print_plitemname (pbp-file, optnr, nd[opt].itm, nd[opt].clr); 
fprintf (pbp-file, "u"); 


} 
printf (pbp_file,"\n"); 
plconstraints ++; 


This code is used in section 84. 
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47. When backtracking, we need to add a RUP constraint. 
(Write RUP constraint on backtrack for proof logging 47) = 


{ 
int opt; 
if (stage >0) { 
{printf (pbp_file, "rup."); 
for (k =0; k < stage; k++) { 
for (opt = choice |levelstage|k]]; nd[opt].itm > 0; opt——) ; 
optnr = linenr (++ opt); 
fprintf (pbp_file, "147 0"O"d.", optnr); 


forintf (pbp_file, ">=,1;\n"); 
pl_constraints ++; 
} 
} 


This code is used in section 75. 


48. When purging options that have no support, we also need to log a RUP constraint. 
(Write RUP constraint for proof logging when detected unsupported option 48) = 


{ 
int optt, k; 
for (optt = opt; nd[optt].itm > 0; optt——) ; 
fprintf (pbp_file, "rupy1y~o"O"d,", linenr (++ optt)); 
for (k =0; k < rup_stage; k++) { 
for (optt = choice [levelstage [k]]; nd[optt].itm > 0; optt——) ; 
forint (pbp_file, "147 0"O"da,", linenr (+4 optt)); 


fprintf (pbp_file, ">=41u;\n"); 
plconstraints +; 


} 


This code is used in section 51. 


49. When we finish our program, we need to notify our proof logging that we should have derived a 
contradiction. 
(Write contradiction line for proof logging 49) = 

{printf (pbp_file, "cu-1\n"); 


This code is used in section 6. 


50. For better performance of the verifier, it is useful to remove derived constraints from deeper levels of 
our backtrack tree, once they become reduntant for the further verification. For each level of our backtrack 
tree, we log when we enter this level. When we backtrack, we wipe all constraints derived on the lower levels, 
since these have become reduntant. 


(Subroutines 7) += 
void write_proof_log_level (int level) 


{ 
fprintf (pbp_file, "#4"O"d\n", level); 
} 
void wipe_proof_log_level (int level) 
{ 


forintf (pbp_file, "w4"O"d\n", level); 
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51. Maintaining supports. It’s time now to implement some of the mechanisms used for the “virtual 
support array S” described earlier. 

First, let’s see what happens when an option goes away. The value returned is 0 if this was the final option 
for some primary item. Otherwise the option’s trigger list will enqueue fixits, to provide replacements for 
any supports that are no longer valid. 

When opt_out deactivates an option, it sets that option’s “age” to cur_age, which measures our progress 
to a complete solution. Options that are purged early, on the basis of fewer assumptions, are “younger” than 
options that are purged later. 

We'll see later that a trigger list may contain hints about the ages of its entries. Such hints are signalled 
by negative entries. 

The opt_out procedure rearranges the entries of a long trigger list by carrying out a “bucket sort,” which 
puts the youngest remaining entries last. This sorting process uses auxiliary arrays trig_head and trig_tail; 
trig-head is assumed to be zero upon entry and exit. 

When we have just included an option and call empty_the_queue, the proof logging needs to also log the 
included option into its RUP constraint. That is why we add a parameter rup_stage, which at that moment 
is stage + 1 instead of stage. 


#define infinite.age (2 * maz_stage + 2) 


(Subroutines 7) += 
int opt_out(int opt,int act,char *typ,int rup_stage) 
i 
register int i, 77, nn, nnp, p, q, qd, pp, ss, t, optp, cutoff, tmin = infinite_age; 
subroutine_overhead ; 
if (vbose & show_purges) { 
forintf (stderr, "\"O"dy"O"syoption,", cur_age, typ); 
print_option(opt + 1, stderr, 0, 1); 
} 
(Delete opt from the sets of all its unpurified items, possibly returning 0 53); 
0, age(opt) = cur_age; 
if (stremp(typ, "blocking")) { 
(Write RUP constraint for proof logging when detected unsupported option 48); 
i 


for (0,p = trigger(opt), pp =0; p; p= pp) { 
o, optp = info(p),q = link(p); 
0, tt = info(q), pp = link (q); 
if (optp < 0) (If all remaining triggers are known to be inactive, set pp = p and break; otherwise 
discard this hint and continue 67); 
(If optp has been deactivated, set t to its age and goto inactive 68); 
(If % isn’t active, set t = cur_age and goto inactive 69); 
000, info(p) = opt, link(q) = fixit (optp); /* change trigger to fixit */ 
if (=fixit(optp)) { /*x we should enqueue optp */ 
0, link (qrear) = getavail(), info(qrear) = optp, grear = link (qrear); 
0, age(optp) = infintte_age; 
} 
0, fixit (optp) = p; 
continue; 
inactive: if (t <0) (Discard this entry and continue 56); 
if (0, trig_-head|[t] = 0) 0, trig_tail|t] = q; 
00, link (q) = trig_head [t], trig_head [t] = p; /* move trigger to temp list t */ 
if (t < tmin) tmin =¢; 
j 


(Replace trigger (opt) by its unused entries, reordered and hinted 70); 
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totopts —; 
return 1; 


} 


52. (Subroutines 7) += 
int purge_the_option (register int opt,int act,char *typ,int rup_stage) 
{ /* opt isn’t at the left spacer */ 
for (opt--; 0, nd[opt].itm > 0; opt——) ; 
return opt_out (opt, act, typ, rup_stage); 


} 


53. After a secondary item has been purified, we mustn’t mess with its set. Secondary items that lie 
between active and the parameter act are in the process of being purified. 


(Delete opt from the sets of all its unpurified items, possibly returning 0 53) = 
for (nn = opt +1; 0, (4 = nd[nn].itm) > 0; nn++) { 
p = nd[nn\.loc; 


if (p > second A (0, pos(it) > act)) continue; /* ti already purified */ 
0, 88 = size(w#z) — 1; 
if (ss =OAp< second) { /* oops: opt was item ii’s only surviving option */ 


if ((vbose & show_details) \ level < show-choices.max / level > maal — show_choices_gap) { 
fprintf (stderr, "\ycan’ tcover"): 
print_item_name (it, stderr); 
fprintf (stderr, "\n"); 


(Clear the queue and return 0 54); 
} 
0, nnp = set [ii + ss]; 
0, size(ii) = 88; 
oo, set[ti + ss] = nn, set|p] = nnp; 
oo, nd[nn].loc = it + ss,nd[nnp].loc = p; 
updates ++; 


i 


This code is used in section 51. 


54. We can’t complete the current options to a viable set that’s domain consistent. So all of the fixit lists 
remaining in the queue must go back into the trigger lists that triggered them. 


(Clear the queue and return 0 54) = 
while (gfront 4 qrear) { 


0,p = gfront, opt = info(p), qfront = link (p), putavail (p); 
(Change the entries of fixit (opt) back to triggers 55 ); 


} 


return 0; 


This code is used in section 53. 
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55. (Change the entries of fixit(opt) back to triggers 55) = 
{ 
for (0, p = fixit(opt); p; p= pp) { 
00, optp = info(p),q = link(p), info(p) = opt; 
0, pp = link (q); /* info(q) is the same for triggers and fixits */ 
00, link (q) = trigger (optp), trigger (optp) = p; 
} 
o, fixit (opt) = 0; 
} 


This code is used in sections 54, 57, and 74. 


56. An option with negative age will never be used, so we needn’t trigger it. 

(I realized later that, in fact, an inactive option with age 0 will also remain inactive. So I could also 
have discarded a few more entries, and used —c instead of —c — 1 in hints. I’ve decided not to make this 
optimization, for fear of breaking something.) 


(Discard this entry and continue 56) = 
‘ 
putavail (p), putavail (q); 
continue; 


} 


This code is used in section 51. 


30 MAINTAINING SUPPORTS XCCDC-PROOF 857 


57. The queue contains options where we’ve left holes in the support matrix. The fixit lists of those options 
tell us where those holes are. 
(Subroutines 7) += 
int empty_the_queue(int rup_stage) 
i 
register int p, q, pp, dq, 8, 88, W, 77, nn, opt, optp; 
subroutine_overhead ; 
while (gfront 4 qrear) { 
0,p = gfront, opt = info(p), gfront = link(p), putavail (p); 
if (fixit(opt) =0) confusion ("queue"); 
if (leak_checking) list.check (0); 
(If opt is no longer active, revert its fixit list and continue 74); 
(Prepare the mark fields for testing compatibility with opt 27); 
for (0, p = fixit(opt); p; p= pp) { 
0,q = link (p); /* ignore info(p), which is irrelevant for now */ 
0, ti = info(q), pp = link(q); /* t is a primary item, not in opt */ 
for (0,5 = ti, ss = s+ size(ii); s< ss; st++) { 
0, optp = set|s]; 
(If optp is compatible with opt, break 28); 


if (s=ss) { /* opt is inconsistent */ 
if (vbose & show_supports) { 
print_option(opt + 1, stderr, —1,1); 
forintf (stderr ,","); 
print.item_name (ii, stderr); 
fprintf (stderr, "\ynotysupported\n"); 
i 
fixit (opt) = p; 
(Change the entries of fixit(opt) back to triggers 55 ); 
if (sopt_out (opt, active, "purging", rup_stage)) return 0; 
break; /* move to another opt */ 
} else (Record optp as the support for opt and ii 61); 
} 
o, fixit (opt) = 0; 
‘ 


return 1; 


} 
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58. Here’s how we get the ball rolling by making every domain consistent in the first place. 
At the beginning, all mark fields are zero. 
(Subroutines 7) += 
int establish_dc(int rup_stage) 
{ 
register int k, aw, 77, nn, opt, optp, p, d, 9, §, 88; 
cur_age = —1; 
gfront = qrear = getavail(); 
for (opt =0; opt < last_node; 0, opt += nd[opt].loc +1) { 
if (leak_checking) list.check (0); 
(Prepare the mark fields for testing compatibility with opt 27); 
for (k =0; k < osecond; k++) { 
0, = item|k]; 
if (0, mark (i) #4 compatstamp) { /* ti not in opt */ 
for (0,5 = ii, ss = s+ size(ii); s< ss; st+) { 
0, optp = set|s]; 
(If optp is compatible with opt, break 28); 
} 
if (s=ss) { /* opt is inconsistent */ 
if (vbose & show_supports) { 
print_option(opt + 1, stderr, —1, 1); 
forintf (stderr,","); 
print_item_name (ii, stderr); 
forintf (stderr ,"jnotysupported\n"); 


if (sopt_out (opt, active, "purging", rup_stage)) return 0; 


break; /* move to the next opt */ 
} else { 

p = getavail(),q = getavail(); 

0, link (p) = 4; 


0, info(q) = wt; 
(Record optp as the support for opt and ii 61); 
} 
} 
} 
i 


return empty_the_queue(rup_stage ); 


} 


59. (Solve the problem 59) = 
d 
totopts = options; 
if (mestablish_dc(0)) { 
if (vbose & show_choices) fprintf(stderr, ""Inconsistentoptions! \n"); 
goto done; 


} 


(Tidy up the initial trigger lists 60); 
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if (vbose & show_choices) fprintf(stderr, ""Initial,consistency,after."O"1ld_mems.\n", mems); 


(Do a backtrack search, maintaining domain consistency 75); 


} 


This code is used in section 2. 
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60. The purged options that appear in trigger lists are useless baggage. 
( Tidy up the initial trigger lists 60) = 
{ 
register int opt, optp, p, q, PP, 44; 
for (opt =0; opt < last_node; o, opt += nd[opt].loc + 1) 
if (0, age(opt) > 0) { 
for (0,p = trigger (opt),qq =—1; p; p= pp) { 
00, optp = info(p),q = link(p), pp = link (q); 
if (0, age(optp) <0) { 
putavail (p), putavail (q); 
if (qq < 0) 0, trigger (opt) = pp; 
else o, link(qq) = pp; 
} else qq = 9; 
} 
} 
} 


This code is used in section 59. 


61. (Record optp as the support for opt and ii 61) = 
{ 
if (vbose & show_supports) { 
print_option(opt + 1, stderr, —1, 1); 
forintf (stderr ,","); 
print_item_name (it, stderr); 
forintf (stderr, "\ysupportedyby.,"); 
print_option(optp + 1, stderr, 0, 1); 


o, info(p) = opt; 
00, link (q) = trigger (optp); 
o, trigger (optp ) = p; 

i 


This code is used in sections 57 and 58. 
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62. A view from the top. Our strategy for generating all exact covers will be to repeatedly choose an 
item that appears to be hardest to cover, namely an item whose set is currently smallest, among all items 
that still need to be covered. And we explore all possibilities via depth-first search, in the following way: 
First we try using the first option in that item’s set; then we explore the consequences of forbidding that 
item. 

The neat part of this algorithm is the way the sets are maintained. Depth-first search means last-in-first- 
out maintenance of data structures; and the sparse-set representations make it particularly easy to undo 
what we’ve done at less-deep levels. 

The basic operation is “covering” each item of a chosen option. Covering means to make an item inactive. 
If it is primary, we remove it from the set of items needing to be covered, and we block all other options 
that contain it. If the item is secondary and still active (not yet purified), we block all options in which it 
has the wrong color. 

The branching discipline that we follow is quite different from what we did in DLX2 or SSXCC, however, 
because we’re now maintaining domain consistency throughout the search. The old way was to choose a 
“best item” p, having say d options, and then to try option 1 of the d possibilities for p, then option 2 of 
those d, ..., option d of those d, before backtracking to the previous level. 

The new way, given consistent domains, starts out the same as before. We choose a best item p;, having 
d, options, and we try its first option. But after returning from that branch, we remove that option and 
restore domain consistency; then we choose a new best item p2, having dz options, and try the first of those. 
Eventually, after trying and removing the first remaining options of p; through p,, we'll reach a point where 
we can’t make the remaining domains both consistent and nonempty. That’s when we back up. 

In this scenario, all of the subproblems for pi, ..., py are trying to extend the same partial solution with 
s choices to a partial solution that has s+ 1 choices. We call this “stage s” of the search. Stage s actually 
involves k different nodes of the (binary) search tree, each of which is on its own “level.” (The level is 
the distance from the root; the stage is the number of options that have been chosen in the current partial 
solution.) 

We might think of the search as a tree that makes a k-way branch at stage s, instead of as a tree that 
makes binary branches at each level. Such an interpretation is equivalent to the “natural correspondence” 
between ordinary trees and binary trees, discussed in TAOCP Section 2.3.2. 
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63. As search proceeds, the current subproblem gets easier and easier as the number of active items and 
options gets smaller and smaller. Let J, be the set of all items that are active when s options c), ..., Cs 
have been chosen to be in the partial solution-so-far. Thus Jo is the set of all items initially given; and I, for 
s > 0 is obtained from J,_; by removing the primary items and the previously unpurified secondary items 
of c,. We denote the primary items of IJ, by P,; these are the primary items not in cj, ..., Cs. 

Let O_, be the set of all options actually given. Just before entering stage 0, we reduce O_, to O}"*, 
the largest subset of O_, that is domain consistent, by purging options that have no support. In general, 
stage s begins with a domain-consistent set of options Oi™*, which is the largest such set that’s compatible 
with c1, ..., cs. Later on in stage s we usually work with a smaller set of active options O;, which is the 
largest domain-consistent set that’s contained in Oi after we’ve removed the options whose consequences 
as potential choices were previously examined in this stage. 

If every item in P, still belongs to at least one option of O,, we’re ready to make a new c,; from among 
those remaining options. We get On from O, by choosing cs; and blocking every option incompatible 
with it, and then by purging options that aren’t domain-consistent. 

Thus when we’re in stage s, there’s a sequence of sets of options 


O42 O05 SOP" S Oj 4 50e" 2.03, 
all of which are domain consistent except possibly O_,. Notice that 
if o € OM and p € o then p € P. 


And there’s good news: The support array S/o, p] follows the nested structure of our search in a useful 
way. Recall that Slo, p| = # if p € 0; otherwise S/o, p] = 0’, where p € o’ and o’ is compatible with o. 

This array is defined for all options o € Oi, and for all primary items p € Py. However, when we enter 
stage s, we’re interested only in the much smaller subarray that contains supports when o € O'™* and p € P,. 
And when we’re transitioning from stage s to stage s + 1, we care only about a still-smaller subarray, for 
o € O, and p € P,. In particular, domain consistency implies that we have 


if o € O@* and p ¢ o and p € P, then Slo, p] € 0"; 
ifo € O, and p g o and pé€ P, then S/o, p] € Os. 


64. Eventually a choice will fail, of course. Backtracking becomes necessary in two distinct ways: (1) If 
we’ve settled on a new c, among the options of O,, but we’re unable to reduce the remaining compatible 
options to a domain-consistent Oe without emptying some domain, we “backtrack in stage s” and reject 
that choice. (Thus, we stay in stage s but move to a new level; the active items remain the same.) (2) If 
we’ve finished exploring a choice from O, and are unable to reduce the other options to a smaller domain- 
consistent O,, we “backtrack to stage s — 1” and reject c,_1. (Thus, we resume where we left off at the 
previous stage’s deepest level; the active items revert back from P, to the larger set P,_1.) 

I wish I could say that it was easy for me to discover the programming logic just described. I guess it was 
my baptism into what researchers have called “fine-grained” versus “coarse-grained” algorithms. 

Notice that when we backtrack, we need not change the S array in any way. A support is always a support. 
Thus there’s no point in trying to undo any of the changes we’ve made to the current support structure. 
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65. The triggering. Suppose there are 1000 options and 100 items. Then the S array has 100,000 
entries, most of which are supports (that is, not #). Every support is an entry in a trigger list; hence the 
trigger lists are necessarily long. The task of maintaining domain consistency might therefore seem hopelessly 
inefficient. 

On the other hand, after we’ve made some choices, there may be only 100 options left, and perhaps 30 
items not yet covered. Then at most 3000 supports are relevant, and most of the information in trigger lists 
is of no interest to us. An efficient scheme might therefore still be possible, if we can figure out a way to 
avoid looking at useless triggers. 

Ideally we’d like options from O, to appear at the top of each trigger stack, with options from Oi just 
below them, and with O,-1, O'™!, ..., Oo, OjM* furthest down. The pairs (0, p) of interest would then 
appear only near the top. 

Unfortunately such an arrangement cannot be guaranteed. Indeed, that’s obvious: The trigger-list entries 
occur in essentially arbitrary order when we first form Oi. If they happen to be supports that work for 
every subsequent stage, no changes to the trigger lists will be needed, and we won’t even want to look at 
those lists. 

We can, however, come sort of close to an ideal arrangement, by exploiting the fact that every option not 
in the current O, has been deactivated at least once. We look at trigger (o) only after o has become inactive; 
and at that time we can reorder its entries. 

Therefore this program inserts markers into the trigger lists, saying that “all further entries of this list are 
young” (meaning deactivated early, hence uninteresting until we’ve backtracked to an early stage). Every 
such marker is accompanied by a time stamp, so that we can recognize later when its message is no longer true. 


66. When deactivating an option from O'* that won’t be in O,, the “current age” cur_age is 2s. And 
when deactivating an option from O, that won’t be in OM it is 2s + 1. 

Thus an inactive option is in O™* if and only if its age is > 2s; and it’s in O, if and only if its age is 
>2s+1. 

Incidentally, I’ve tried to avoid making bad puns based on cur_age versus courage, or age versus stage. 


67. A negative entry optp = —c in a trigger list is a hint that all future entries will have age less than c. 
The search tree may have changed since this hint was put into the list; so we must look at the relevant stage 
stamp, to ensure that the hint is still valid. 

Suppose o has age 2s. Then o is in O'™* but not in O,. As computation proceeds, without backtracking 
to stage s — 1, the set O, might get smaller and smaller, but o will still not be in O,. Therefore a trigger 
hint saying that o is inactive will be valid until stagestamp[s] changes. (More precisely: If we backtrack to 
stage s — 1, stagestamp[s] won’t change until we progress again to stage s; before that time, we won’t be 
looking at the hint.) 

Suppose o has age 2s+1. Then o is in O, but not in Oa As computation proceeds, without backtracking 
in or to stage s, the set om won’t change. Therefore a trigger hint saying that o is inactive will be valid 
until stagestamp[s + 1] changes. 

That’s why the following code says ‘(cutoff + 1) >> 1’ when selecting the relevant stage stamp. 


(If all remaining triggers are known to be inactive, set pp = p and break; otherwise discard this hint and 
continue 67) = 
{ 


cutoff = —optp — 1; 
if (cutoff < cur_age A (0, ii = stagestamp|(cutoff +1) > 1])) { 
PP =P, 
break; 
} 
putavail (p), putavail (q); /* discard an obsolete hint */ 
continue; /* and ignore it */ 


} 


This code is used in section 51. 
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68. If optp is inactive, it has been purged and its recorded age is cur_age or less. Thus we can conclude 
that optp is active whenever age(optp) > cur_age. 

In general, of course, that age test won’t be conclusive and a slightly more expensive test needs to be made 
by looking further into the data structures. Option optp is active if and only if it appears in the current set 
of its first item. (This is where we use the fact that the first item of optp is primary.) 

(If optp has been deactivated, set t to its age and goto inactive 68) = 

0, t = age(optp); 

if (t < cur_age) { 

0,7) = ndloptp + 1].itm; /* jj is optp’s first item */ 
if (0, nd[optp + 1].loc > jj + size(jj)) goto inactive; 
} /* branch if optp was removed from jj’s set */ 


This code is used in section 51. 


69. When the trigger list for opt refers to an item 7, that item is in opt. Suppose 7 is currently inactive; 
then we wouldn’t be purging opt unless ii has just become inactive (and we’re calling opt_out from within 
include_option). 


(If # isn’t active, set t = cur_age and goto inactive 69) = 
if (0, pos(ii) > active) { 
if (pos(i#) > act) confusion("active"); 
t = cur_age; 
goto inactive; 


} 


This code is used in section 51. 


70. When we get here, pp is either zero or the cell where we found cutoff. In the latter case, pp = p and 
link (p) = q; thus the cutoff hint is in p and q. 
All of the unused trigger entries have been redirected to the trig_head lists, sorted by their age. 


(Replace trigger (opt) by its unused entries, reordered and hinted 70) = 
if (pp =0) cutoff = —-1; 
if (tmin < cutoff) { 
if (tmin < cutoff) confusion ("trig"); 
0, pp = link (q), putavail (p), putavail (q); /* avoid double hint */ 
} 
for (t = tmin; t < cur_age; t++) 
if (0, trig_head[t]) { 
00, link (trig_tail[t]) = pp; 
0, p = getavail(),q = getavail(), link (p) = 4; /* make new hint «/ 
0, info(p) = -t — 1, 
00, info(q) = stagestamp|(t + 1) > 1], link(q) = trig_head [t]; 
o, trig_head[t] = 0; 
PP =P; 


if (trig_head[cur_age]) { 
00, link (trig_tail[cur_age]) = pp; /* give no hint for inactive options of the current age «/ 
0, pp = trig-head[cur_age], trig_-head |cur_age] = 0; 

} 

0, trigger (opt) = pp; 


This code is used in section 51. 
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71. (Print a trigger hint 71) = 
4 
fprintf (stderr, "cutoff ,forjage,"O"d", —info(p) — 1); 
if (info(q) 4 stagestamp |(—info(p)) >> 1]) fprintf (stderr, "\,Cobsolete)"); 


This code is used in section 21. 


72. At this point we want curstamp to have a value that’s larger than anything found in a trigger list hint. 
Moreover, the values of stagestamp [0], ..., stagestamp|stage —1] should all be distinct and less than curstamp, 
because they might be used in future hints. 

We may not be able to satisfy those conditions when badstamp is a small positive constant! But we will 
have checked out the following code at least once before failing. 


(Bump curstamp 72) = 
if (++ biggeststamp = badstamp) { 
if (badstamp > 0A stage > badstamp) { 
forintf (stderr, "Timestamp overflow,,(badstamp="O"d) !\n", badstamp); 
exit (—11); 
i 
(Remove all hints from all trigger lists 73); 
for (k = 0; k < stage; k++) 0, stagestamp|k] = k; 
biggeststamp = k; 
} 
curstamp = biggeststamp; 


This code is used in section 77. 


73. Therefore, when curstamp “wraps around,” we must abandon all of the hints that were to be validated 
by obsolete timestamps. 
(Remove all hints from all trigger lists 73) = 
for (k =0; k < last_node; k += nd[k].loc +1) { 
for (p = trigger(k); p; 0,p = link(p)) 
if (0, info(p) <0) { 
0,q = link(p),r = link (q); /* we know that link(q) #0 */ 
00, info(p) = info(r), link (p) = link (r); 
putavail (q), putavail (r); 
} 
} 


This code is used in section 72. 


74. When opt was put into the queue, we made its age infinite. So it will have been purged in the meantime 
if and only if its age is now cur_age. 
(If opt is no longer active, revert its fixit list and continue 74) = 
if (0, age(opt) 4 infinite_age) { 
(Change the entries of fixit (opt) back to triggers 55 ); 
continue; 


} 


This code is used in section 57. 
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75. The dancing. 


(Do a backtrack search, maintaining domain consistency 75) = 
level = stage = —1; 
newstage: (Increase stage 77); 
newlevel: nodes ++; 
(Increase level 78); 
write_proof_log_level (level ); 
if (vbose & show_profile) profile |stage]++; 
if (sanity_checking) sanity(); 
if (leak_checking) list.check (0); 
(Do special things if enough mems have accumulated 79); 
if (stage < groundstage) (Read and act on an option from zcutoff_file 94); 
tmems = mems; 
if (vbose & show_option_counts) fprintf(stderr, 
"Level "O"d, ystage,,"O"d,,"O"d options, "O"d items\n", level, stage, totopts, active); 
(Set best_itm to the best item for branching and t¢ to its size 83); 
bmems += mems — tmems; 
if (stage = xcutoff) (Output a partial solution and goto backup 91); 
if (t = infty) (Visit a solution and goto backup 84); 
00, choice [level] = cur_choice = set[best_itm]; 
got_choice: o, deg|level] = t; 
if (t = 1) 0, saved|stage] = saveptr; 
else (Save the currently active sizes 85); 
cur_age = stage + stage + 1; 
tmems = mems; 
if (sinclude_option(cur-choice)) { 
nmems += mems — tmems; 
goto tryagain; 


if (sempty_the_queue(stage +1)) { 
nmems += mems — tmems; 
goto tryagain; 
} 
goto newstage; 
tryagain: if ({ =1) goto backup; 
if (vbose & show_choices) fprintf(stderr, "Backtracking,,in,stage,"O"d\n", stage); 
goto purgeit; 
backup: write_proof_log_level (levelstage [stage — 1]); 
(Write RUP constraint on backtrack for proof logging 47); 
if (deg [level] 4 1 A levelstage[stage — 1] #0) { 
wipe_proof_log_level (level ); 
} 
if (—stage < groundstage) goto done; 
if (vbose & show_choices) fprintf(stderr, "Backtracking,to,stage,,"O"d\n", stage); 
o, level = levelstage [stage]; 
if (0, deg[level] =1) goto backup; 
purgeit: if (ubose & show_option_counts) totopts = levelopts [level]; 
(Restore the currently active sizes 86 ); 
new_age: cur_age = stage + stage; 
tmems = mems; 
if (A(0, purge_the_option (choice [level], active, "removing", stage))) { 
pmems += mems — tmems; 
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goto backup; 


} 

if (sempty_the_queue(stage)) { 
pmems += mems — tmems; 
goto backup; 


} 


goto newlevel; 


This code is used in section 59. 


76. We save the sizes of active items on savestack, whose entries have two fields | and r, for an item and 
its size. This stack makes it easy to undo all deletions, by simply restoring the former sizes. 


(Global variables 4) += 


int stage; /* number of choices in current partial solution «/ 

int level; /* current depth in the search tree (which is binary) */ 

int cur_age; /* current stage or stage + 1/2 (times 2) «/ 

int choice|maz_level]; /* the option and item chosen on each level */ 

int deg|maz_level]; /* the number of options that item had at the time */ 
int saved[maa_stage]; /* size of savestack at each stage */ 

int levelstage|max_stage]; /* the most recent level at each stage */ 

int stagelevel |maz_level]; /* the stage that corresponds to each level «/ 
int levelopts|maz_level]; /* options remaining at each level */ 

int stagestamp|maz_stage]; /* timestamp that’s current at each stage «/ 
ullng profile |maz_stage] = {1}; /* number of search tree nodes on each stage */ 
twoints savestack |savesize]; 

int saveptr; /* current size of savestack */ 

int trig_head [infinite_age], trig_tail[infinite_age]; /* in opt_out */ 


77. (Increase stage 77) = 
if (++stage > mars) { 
if (stage > maa_stage) { 
fprintf (stderr, "Too many,stages !\n"); 
exit (—40); 
} 


mazs = stage; 
} 
(Bump curstamp 72); 
o, stagestamp|stage] = curstamp; 


This code is used in section 75. 


78. (Increase level 78) = 
if (++level > maal) { 
if (level > maz_level) { 
forintf (stderr, "Too,many,levels!\n"); 
exit (—4); 
} 


mazl = level; 


} 


oo, stagelevel [level] = stage, levelstage|stage| = level; 
if (vbose & show_option_counts) levelopts [level] = totopts; 


This code is used in section 75. 
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79. (Do special things if enough mems have accumulated 79) = 
if (delta A (mems > thresh)) { 
thresh += delta; 
if (vbose & show_full_state) print_state( ); 
else print_progress(); 
} 
if (mems > timeout) { 
fprintf (stderr, "TIMEOUT! \n"); goto done; 


} 


This code is used in section 75. 


80. This is where we extend the partial solution. 
Notice a tricky point: We must go through the sets from right to left, because the options we block move 
right as they leave the set. 
(Subroutines 7) += 
int include_option(int opt) 
{ 
register int c, cc, k, p, q, pp, 8, 88, optp; 
subroutine_overhead ; 
if (vbose & show_choices) { 
forintf (stderr, "S"O"d:", stage); 
print_option (opt, stderr, 1,0); 
} 
for (opt—-; 0, nd[opt].itm > 0; opt——) ; /* move back to the spacer */ 
(Inactivate all items of opt, and record their colors 81 ); 
for (k = active; k < oactive; k++) { 
00,8 = item|k], ss = s + size(s) — 1; 
if (s > second A (0,c = match(s))) { /* we must purify s */ 
for (; ss >s; ss——) { 
0, optp = set|ss]; 
if ((0, nd[optp].clr 4 c) A apurge_the_option (optp, oactive, "blocking",0)) return 0; 


} else 
for (; ss > s; ss——) { 
0, optp = set|ss] — 1; 
while (0, nd[optp].itm > 0) optp—-; /* move to the spacer */ 
if (optp 4 opt \ sopt_out(optp, oactive, "blocking",0)) return 0; 
} 
} 
(Make opt itself inactive 82); 
return 1; 


} 
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81. An item becomes inactive when it becomes part of the solution-so-far (hence it leaves the problem- 
that-remains). Active primary items are those that haven’t yet been covered. Active secondary items are 
those that haven’t yet been purified. 

The active items are the first active entries of the item list. At one time I thought it would be wise to 
keep primary and secondary items separate, using a sparse-set discipline independently on each sector. But 
I found that the time spent in maintaining and searching the active list was negligible in comparison with 
the overall running time; so I’ve kept the implementation simple. 

At this point in the computation, an item of opt will be inactive if and only if it is secondary and purified, 
because we are including opt in the partial solution. 


(Inactivate all items of opt, and record their colors 81) = 
p = oactive = active; 
for (q = opt +1; 0,(c = nd[qj.itm) > 0; q+) { 
0, pp = pos(c); 
if (pp <p) { /* cis active */ 
0, cc = item|{——p}; 
00, item|p] = c, item[pp] = cc; 
if (c > second) 00, match(c) = nd{q].clr; 
00, pos(cc) = pp, pos(c) = p; 
updates ++; 
} 
} 


active = p; 


This code is used in section 80. 


82. This program differs from SSXCC in one significant way: It makes option opt inactive. In particular, it 
makes all of opt’s items have size 0, except for unpurified secondaries. (Thus, we essentially say that newly 
assigned variables—the inactivated primary items—should have empty domains when they leave the current 
subproblem, while SSXCC left them with domains of size 1.) 

It would be a mistake to call opt_out (opt, oactive), of course, because that procedure doesn’t want any 
primary items to become optionless. On the contrary, we have precisely the opposite goal: We celebrate the 
fact that all of the primaries in opt have become covered. 

We don’t have to change trigger (opt), because no active options involve inactive primary items. 


(Make opt itself inactive 82) = 

for (k = active; k < oactive; k++) { 
0, 8 = item|k]; 
if (s > second) continue; 
if (size(s) £1) confusion("include"); 
0, size(s) = 0; 

} 

if (vbose & show_purges) { 
forintf (stderr, "\"O"d choosing, option,.", cur_age); 
print_option (opt + 1, stderr, 0,1); 


0, age(opt) = cur_age; 
totopts —— ; 


This code is used in section 80. 
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83. The “best item” is considered to be an item that minimizes the number of remaining choices. If there 
are several candidates, we choose the first one that we encounter. 

Each primary item should have at least one valid choice, because of domain consistency. 
#-define infty *7£ffffff 


(Set best_itm to the best item for branching and ¢ to its size 83) = 
t = infty; 
if ((vbose & show-_details) \ level < show_choices_max / level > maal — show_choices_gap ) 
fprintf (stderr, "Stage,"O"d,", stage); 
for (k =0; t>1Ak < active; k++) 
if (0, item[k] < second) { 
0, 8 = size(item|k]); 
if ((vbose & show_details) A level < show choices.max / level > maal — show-choices.gap) { 
print_item_name (item|[k], stderr); 


fprintf (stderr, "("O"d)", s); 


ws <t) { 
f (s =0) fprintf (stderr, "I’m confused. \n"); /* hide missed this */ 
if (s < t) best_itm = item|k],t = s; 

Sits if (item|k] < best_itm) best_itm = item |k]; 


if ((vbose & show_details) A level < show choices.max / level > maal — show-choices_gap) { 
if (t = infty) fprintf (stderr ,"jsolution\n"); 
else { 
fprintf (stderr, "\joranching,,on"); 
print_item_name (best_itm, stderr ); 


forintf (stderr, "("O"d) \n",t); 


if (t > mardeg At < infty) mardeg =t; 
if (shape_file) { 
if (¢ = infty) fprintf (shape_file,"sol\n"); 
else { 
forintf (shape_file,""O"d", t); 
print_item_name (best_itm, shape_file); 
forintf (shape_file, "\n"); 


fflush(shape_file); 


This code is used in section 75. 
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84. (Visit a solution and goto backup 84) = 
{ 

count ++; 

( Write a solution for proof logging 46); 

if (spacing A (count mod spacing =0)) { 
printf (""O"11d:\n", count); 
for (k =0; k < stage; k++) print_option (choice [levelstage [k]], stdout, 0, 0); 
fflush (stdout); 


if (count > maxcount) goto done; 
goto backup; 


} 


This code is used in section 75. 


85. (Save the currently active sizes 85) = 
{ 
o, saved [stage] = saveptr; 
if (saveptr + active > maxsaveptr) { 
if (saveptr + active > savesize) { 
fprintf (stderr, "Stack,overf low, (savesize="O"d) !\n", savesize ); 
exit (—5); 
} 
maxsaveptr = saveptr + active; 
} 
for (p = 0; p < active; p++) 
000, savestack|saveptr + p].l = item|[p], savestack|saveptr + p].r = size (item[p]); 
saveptr += active; 


} 


This code is used in section 75. 


86. (Restore the currently active sizes 86) = 
o, active = saveptr — saved |stage]; 
saveptr = saved |stage]; 
for (p = 0; p < active; pt+) 00, size(savestack|saveptr + p].l) = savestack|saveptr + p].r; 


This code is used in section 75. 


87. (Subroutines 7) += 
void print_savestack (int start, int stop) 


{ 
register int k; 
for (k = start; k < stop; k++) { 
print_item_name (savestack [k].1, stderr); 
forintf (stderr, "("O"d) ,"O"d\n", savestack [k].1, savestack [k].r); 
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88. (Subroutines 7) += 
void print_state (void) 


{ 
register int I, s; 
fprintf (stderr, "Currentystate,(level,"O"d, ystage,"O"d) :\n", level, stage); 
for (1 =0; | < level; I++) { 
if (levelstage|stagelevel [l]] A 1) fprintf(stderr, "~"); 
print_option (choice [I], stderr, —1, 1); 
fprintf (stderr, "\Cofy"O"d)", deg [l]); 
if (vbose & show_option_counts) fprintf (stderr ,",4"O"d opts\n", levelopts[l]); 
else fprintf(stderr ,"\n"); 
if (1 > show_levels.max) { 
forintf (stderr, "\y...\n"); 
break; 
} 


forintf (stderr, "\y"O"11ldjsolutions,,"O"1ld.mems, ,,and,max,level,"O"dso yfar.\n", count, 
mems, mazl); 
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89. During a long run, it’s helpful to have some way to measure progress. The following routine prints a 
string that indicates roughly where we are in the search tree. The string consists of node degrees, preceded 
by ‘~’ if the node wasn’t the current node in its stage (that is, if the node represents an option that has 
already been fully explored — “we’ve been there done that”). 

Following that string, a fractional estimate of total progress is computed, based on the naive assumption 
that the search tree has a uniform branching structure. If the tree consists of a single node, this estimate is .5. 
Otherwise, if the first choice is the kth choice in stage 0 and has degree d, the estimate is (k —1)/(d+k-1) 
plus 1/(d+k-—1) times the recursively evaluated estimate for the kth subtree. (This estimate might obviously 
be very misleading, in some cases, but at least it tends to grow monotonically.) 

Fine point: If we’ve just backtracked within stage stage, the string of node degrees with end with a ‘~ 
entry, and we haven’t yet made any choice in the current stage. The test ‘l = level — 1’ below uses the fact 
that levelstage[stage] = level to adjust the fractional estimate appropriately for the partial progress in the 
current stage. 

(Subroutines 7) += 
void print_progress (void) 


{ 


2 


register int J, Il, k, d, c, p, ds =0; 
register double f/f, jd; 
fprintf (stderr, "yafter,"O"1lldmems:"O"1lldjsols,", mems, count); 
if (stage < groundstage) 
forintf (stderr, "\yinitializing,atstage,"O"d,level,"O"d\n", stage, level); 
else { 
for (f = 0.0, fd = 1.0,1 = (groundstage ? levelstage|groundstage — 1] + 1:0); | < level; I++) { 
if (I < show_levels_maxz + levelstage |groundstage]) 
forintf (stderr, "\"O"s"O"d", levelstage [stagelevel [l]] =1?"" :"~", deg[l]); 
if (levelstage [stagelevel [l]] =1V1 = level —1) { /* see remark above */ 
for (k = 1,d= deg[l], ll =1—1; Il > 0A stagelevel [ll] = stagelevel|l]; kt++,d++,ll—) ; 
fd x= d, f += (k—-1)/fa; /* choice | is treated like k of d «/ 


if (1 > show_levels_maxz + levelstage|groundstage] \ ds) ds = 1, fprintf(stderr,"..."); 
} 
forintf (stderr, ""O".5£\n", f + 0.5/fd); 


} 


90. (Print the profile 90) = 
{ 
fprintf (stderr, "Profile:\n"); 
for (k = 0; k < mazs; k++) fprintf(stderr ,""O"3d:\4"O"11d\n", k, profile [k]); 
} 


This code is used in section 6. 
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91. I’m experimenting with a mechanism by which partial solutions of a large problem can be saved to 
temporary files and computed separately — for example, by a cluster of computers working in parallel. Each 
partial solution can be completed to full solutions when this program is run with one of the files output here, 
using X( filename ) on the command line. 


(Output a partial solution and goto backup 91) = 


(Open a new zcutoff_file 92); 

fprintf (acutoff_file, "Resumeatystage,"O"d\n", stage); 

for (k =0; k < level; k++) { 
for (0,7 = choice|k]; 0, nd[j — 1].ttm > 0; j--) ; 
putc(levelstage [stagelevel[k]] Ak ? ?-? : 2+’, xcutoff_file); 
print_option(j, xcutoff_file, 0, 1); 

} 

fclose (acutoff_file ); 

xcount ++; 

goto backup; 


} 


This code is used in section 75. 


92. #define part_file_prefix "“/tmp/part" /* should be at most 10 or so characters */ 
#define part_filename_size 20 


(Open a new xcutoff_file 92) = 
k = sprintf (xcutoff_name, part_file_prefixO"da" , zcount); 
xcutoff_file = fopen(a«cutoff_name, "w"); 
if (s2cutoff_file) { 
fprintf (stderr, "Sorry, ,Ican’ tyopen file,‘"O"s’ for writing! \n", ccutoff_name); 
exit (—1); 
} 


This code is used in section 91. 


93. (Open scutoff_file for reading, and break 93) = 
strnepy (acutoff_name, argu [j] + 1, part_filename_size — 1); 
xcutoff_file = fopen(axcutoff_name,"xr"); 
if (nxcutoff_file ) 
forintf (stderr, "Sorry, ,Ijcan’ tyopen file,‘"O"s’ for reading! \n", ccutoff_name); 
if (fgets (buf , bufsize, xcutoff_file)) { 
if (strncmp (buf, "Resumejatstage,,",16) = 0) { 
for (groundstage = 0,1 = 16; 0, buf [i] > °0’ A buf [é] < 797; i+4) 
groundstage = 10 * groundstage + buf [i] — ’0’; 
if (vbose & show_basics) fprintf(stderr, "Resuming atystage,,"O"d\n", groundstage ); 
} 


} 
break; 


This code is used in section 5. 
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94. (Read and act on an option from zcutoff_file 94) = 

if (fgets (buf , bufsize, xcutoff_file)) confusion ("resuming"); 

0, choice [level] = cur_choice = read_option (); 

if (cur_choice <0) { 
forintf (stderr, "Misformatted option in file, ‘"O"s’:\n", ccutoff_name); 
forintf (stderr,""O"s", buf ); 
exit (—1); 

} 


t=1; 
if (0, buf [0] = ’+’) goto got_choice; 
goto neu_age; 


} 


This code is used in section 75. 


95. (Report the number of partial solutions output 95) = 
forintf (stderr, "Partial.solutions,saved,on,," part_file_prefiz"0.."O"s.\n", xccutoff_name); 


This code is used in section 6. 


96. (Global variables 4) += 
char xcutoff_name |part_file.name_size]; 
FILE x«cutoff_file ; 
int groundstage; /* the stage where calculation begins or resumes */ 
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97. Index. 
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f: 89. 
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fflush: 
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got.choice: 75, 94. 
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44, 46, 47, 48, 52, 53, 68, 80, 81, 91. 

Cay 

ii: 22, 28, 51, 57, 58, 68. 
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last.itm: 10, 29, 30, 32, 34. 

last.node: 6, 10, 13, 16, 22, 25, 31, 32, 33, 36, 
38, 42, 58, 60, 73. 

leak.checking: 6, 25, 57, 58, 75. 
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levelopts: 75, 76, 78, 88. 

levelstage: 46, 47, 48, 75, 76, 78, 84, 88, 89, 91. 
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58, 60, 61, 70, 73. 
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il: 89. 
Iname: 9, 12, 14, 29, 30, 32, 35, 41. 


loc: 9, 13, 16, 22, 25, 31, 32, 36, 53, 58, 60, 68, 73 
Ir: 11, 12, 14, 29, 30, 31, 32, 41. 


main: 2. 

mark: 9, 26, 27, 28, 35, 58. 
match: 9, 27, 28, 80, 81. 
max.cols: 2, 10, 29, 34. 
maz_level: 2, 76, 78. 
max.nodes: 2, 10, 31, 32. 
mazstage: 2, 51, 76, 77. 
maxcount: 4, 5, 84. 
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maal: 4, 6, 53, 78, 83, 88. 

mars: A, 6, 77, 90. 

maxsaveptr: 4, 6, 85. 

2, 3, 4, 6, 7, 25, 59, 75, 79, 88, 89. 
2, 84. 
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mod: 
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namebuf : 

nb: 35. 

nd: 9, 10, 13, 14, 16, 19, 22, 25, 27, 28, 31, 
32, 33, 36, 42, 44, 46, 47, 48, 52, 53, 58, 60, 
68, 73, 80, 81, 91. 

new_age: 75, 94. 

newlevel: 75. 

newstage: 7d. 

nmems: A, 6, 75. 

mn: 27, 28, 51, 53, 57, 58. 

nnp: 51, 53. 

node: 6, 9, 10. 

node_struct: 9. 

nodes: 4, 6, 7, 75. 

O: 2. 

o: 2. 

oactive: 10, 80, 81, 82. 

oo: 2, 14, 26, 29, 32, 33, 34, 35, 51, 53, 55, 60, 
61, 70, 73, 75, 78, 80, 81, 86. 

ooo: 2, 51, 85. 

opb_file: 4, 5, 6, 40, 42, 43, 44. 

opbname: 4, 5. 

opt: 13, 14, 19, 21, 22, 23, 25, 27, 28, 43, 46, 
A7, 48, 51, 52, 53, 54, 55, 57, 58, 60, 61, 69, 
70, 74, 80, 81, 82. 


11, 12, 14, 29, 30, 31, 32, 41, 42. 


opt_out: 51, 52, 57, 58, 69, 76, 80, 82. 
options: A, 31, 38, 59. 

optnr: 2, 41, 42, 46, 47. 

optp: 22, 28, 51, 55, 57, 58, 60, 61, 67, 68, 80. 
optt: A8. 

opti: AA. 

opti: AA4. 

opt2: AA. 

osecond: 10, 34, 35, 38, 58. 


p: 2,13, 15, 20, 21, 23, 24, 25, 51, 57, 58, 60, 80, 89. 

panic: 29, 30, 31, 32. 

part_filename_size: 92, 93, 96. 

part_file_prefiz: 92, 95. 

pbp-_file: 4, 5, 6, 45, 46, 47, 48, 49, 50. 

pbp_name: 4, 5. 

plLconstraints: 4, 40, 42, 43, 44, 45, 46, 47, 48. 

pmems: A, 6, 75. 

pool: 2, 18, 19, 20, 25. 

poolptr: 6, 19, 20, 25. 

poolsize: 2, 19, 20. 

pos: 9, 15, 16, 32, 33, 35, 53, 69, 81. 

pp: 2, 31, 32, 51, 55, 57, 60, 67, 70, 80, 81. 

primextra: 9, 15, 34. 

print_fixit: 23. 

printitem_name: 12, 13, 15, 16, 21, 23, 37, 53, 
57, 58, 61, 83, 87. 

printitm: 15. 
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print_option: 13, 14, 21, 23, 24, 51, 57, 58, 61, 
80, 82, 84, 88, 91. 

print_plitemname: 41, 42, 43, 44, 46. 

print_progress: 79, 89. 

print_queue: 24. 

print_savestack: 87. 

print_state: 79, 88. 

printtrigger: 21, 22. 

print_triggers: 22. 

printf: 84. 

profile: 6, 75, 76, 90. 

propt: 13. 

prow: 13, 15. 

purge_the_option: 

purgeit: 75. 

putavail: 20, 54, 56, 57, 60, 67, 70, 73. 

putc: 91. 

q: 2, 13, 14, 16, 21, 23, 51, 57, 58, 60, 80. 

qfront: 19, 24, 25, 54, 57, 58. 

qq: 16, 28, 51, 57, 58, 60. 


52, 75, 80. 


grear: 19, 24, 25, 51, 54, 57, 58. 

r: 2, ll, 16. 

read_option: 14, 94. 

rname: 9, 12, 14, 29, 30, 32, 35, 41. 
rup_stage: 48, 51, 52, 57, 58. 

s: 2, 57, 58, 80, 88. 

sanity: 16, 75. 

sanity.checking: 6, 16, 75. 

saved: 75, 76, 85, 86. 

saveptr: 75, 76, 85, 86. 

savesize: 2, 76, 85. 

savestack: 2, 4, 76, 85, 86, 87. 
second: 10, 13, 15, 27, 28, 29, 31, 32, 34, 35, 39, 


41, 43, 44, 53, 80, 81, 82, 83. 
secondextra: 9, 34. 
set: 9, 10, 11, 15, 16, 19, 29, 35, 36, 43, 44, 
53, 57, 58, 75, 80. 
setlength: 6, 10, 15, 34. 
shape_file: 4, 5, 6, 83. 
shape.name: 4, 5. 


show_basics: 2, 3, 4, 6, 93. 
show_choices: 3, 4, 37, 59, 75, 80. 
show_choices.gap: 4, 5, 53, 83. 
show_choices_maz: 4, 5, 53, 83. 
show_details: 3, 4, 53, 83. 
show-_full_state: 3, 79. 
show_levels_mazr: 4, 5, 88, 89. 


show-maz_deg: 3, 6. 
show_mstats: 3, 6. 
show_option_counts: 3, 75, 78, 88. 
show-profile: 3, 6, 75. 
show_purges: 3, 51, 82. 
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show_supports: 3, 57, 58, 61. vbose: 2, 3, 4, 5, 6, 31, 37, 51, 53, 57, 58, 59, 61, 
show_tots: 2, 3. 75, 78, 79, 80, 82, 83, 88, 93. 
show.warnings: 3, 4, 31. vet.and_set: 25. 
showid: 18. wipe_proof_log_level: 50, 75. 
showpos: 13. write_proof_log_level: 50, 75. 
signbit: 25. z: 13, 14, 16. 
size: 9, 13, 15, 16, 22, 32, 33, 34, 35, 39, 43, 44, zcount: 4, 6, 91, 92. 

53, 57, 58, 68, 80, 82, 83, 85, 86. zcutoff: 4, 5, 75. 
spacing: 4, 5, 84. xcutoff_file: 91, 92, 93, 94, 96. 
sprintf: 92. xcutoff_name: 92, 93, 94, 95, 96. 
8S: 43, 44, 51, 53, 57, 58, 80. xtra: 9, 19. 
sscanf: 5. 


stage: 46, 47, 51, 72, 75, 76, 77, 78, 80, 83, 84, 
85, 86, 88, 89, 91. 

stagelevel: 6, 76, 78, 88, 89, 91. 

stagestamp: 67, 70, 71, 72, 76, 77. 

start: 87. 

stderr: 2, 3, 5, 6, 7, 18, 15, 16, 20, 21, 23, 24, 
25, 29, 31, 37, 38, 39, 51, 53, 57, 58, 59, 61, 
71, 72, 75, 77, 78, 79, 80, 82, 83, 85, 87, 88, 
89, 90, 92, 93, 94, 95. 


stdin: 29, 31. 

stdout: 84. 

stop: 87. 

str: 11, 12, 14, 29, 31, 41, 42. 
streat: 5. 

stremp: 51. 

strepy: 5. 


stream: 12, 13, 41. 
stringbuf: 11. 

strlen: 29, 31. 

strncmp: 93. 

strncpy: 93. 
subroutine_overhead: 2, 51, 57, 80. 
$1: 4a. 

62: 44. 

t: 2, 25, dl. 

thresh: 4, 5, 79. 

temeout: 3, 4, 5, 79. 

tmems: 4, 75. 

tmin: 51, 70. 

totopts: 10, 51, 59, 75, 78, 82. 
trig-head: 51, 70, 76. 

trig-tail: 51, 70, 76. 

trigger: 19, 21, 25, 51, 55, 60, 61, 70, 73, 82. 
tryagain: 75. 

twoints: 6, 11, 19, 76. 

typ: 51, 52. 

uint: 2. 

ullng: 2, 4, 76. 

updates: 4, 6, 53, 81. 
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(Adjust nd 36) Used in section 31. 

(Bump compatstamp 26) Used in section 27. 

(Bump curstamp 72) Used in section 77. 

(Change the entries of fixit (opt) back to triggers 55) Used in sections 54, 57, and 74. 
(Check for duplicate item name 30) Used in section 29. 

(Clear the queue and return 0 54) Used in section 53. 

(Create a node for the item named in buf [p] 32) Used in section 31. 

(Delete opt from the sets of all its unpurified items, possibly returning 0 53) Used in section 51. 
(Discard this entry and continue 56) Used in section 51. 

(Do a backtrack search, maintaining domain consistency 75) Used in section 59. 
(Do special things if enough mems have accumulated 79) Used in section 75. 
(Expand set 35) Used in section 31. 

(Global variables 4, 10, 19, 76, 96) Used in section 2. 

( 
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If all remaining triggers are known to be inactive, set pp = p and break; otherwise discard this hint and 


continue 67) Used in section 51. 
(If ii isn’t active, set t = cur_age and goto inactive 69) Used in section 51. 
(If optp has been deactivated, set t to its age and goto inactive 68) Used in section 51. 
(If optp is compatible with opt, break 28) Used in sections 57 and 58. 
(If opt is no longer active, revert its fixit list and continue 74) Used in section 57. 
(Inactivate all items of opt, and record their colors 81) Used in section 80. 
(Increase level 78 ) Used in section 75. 
(Increase stage 77) Used in section 75. 
(Initialize item 34) Used in section 31. 
(Input the item names 29) Used in section 2. 
(Input the options 31) Used in section 2. 
(Make opt itself inactive 82) Used in section 80. 
(Open a new zcutoff_file 92) Used in section 91. 
(Open «cutoff_file for reading, and break 93) Used in section 5. 
(Output a partial solution and goto backup 91) Used in section 75. 
(Output constraints for options and items 42) Used in section 40. 
(Output constraints to ensure each primary item is covered 43) Used in section 40. 
(Output constraints to ensure secondary items are assigned only one color 44) Used in section 40. 
(Output header of the PBP file for proof logging 45) Used in section 2. 
(Output our problem to the OPB file for proof logging 40) Used in section 2. 
(Prepare the mark fields for testing compatibility with opt 27) Used in sections 57 and 58. 
(Print a trigger hint 71) Used in section 21. 
(Print the profile 90) Used in section 6. 
(Process the command line 5) Used in section 2. 
(Read and act on an option from xcutoff_file 94) Used in section 75. 
(Record optp as the support for opt and ii 61) Used in sections 57 and 58. 
(Remove all hints from all trigger lists 73) Used in section 72. 
(Remove last_node from its item list 33) Used in section 31. 
(Replace trigger (opt) by its unused entries, reordered and hinted 70) Used in section 51. 
(Report an uncoverable item 37) Used in section 2. 
(Report the item totals 39) Used in section 2. 
(Report the number of partial solutions output 95) Used in section 6. 
(Report the successful completion of the input phase 38) Used in section 2. 
(Restore the currently active sizes 86) Used in section 75. 
(Save the currently active sizes 85) Used in section 75. 
(Say adieu 6) Used in section 2. 
(Set best_itm to the best item for branching and t¢ to its size 83) Used in section 75. 
(Solve the problem 59) Used in section 2. 
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Subroutines 7, 12, 13, 14, 15, 16, 20, 21, 22, 23, 24, 25, 41, 50, 51, 52, 57, 58, 80, 87, 88, 89) Used in section 2. 
Tidy up the initial trigger lists 60) Used in section 59. 

Type definitions 9,11) Used in section 2. 

Visit a solution and goto backup 84) Used in section 75. 

Write RUP constraint for proof logging when detected unsupported option 48) Used in section 51. 
Write RUP constraint on backtrack for proof logging 47) Used in section 75. 

Write a solution for proof logging 46) Used in section 84. 

Write contradiction line for proof logging 49) Used in section 6. 
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